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ABSTRACT 

This  report  establishes  a  formal  framework  for  the  kinematic  analysis  of  mechanical  devices. 
The  possible  kinematic  behaviors  of  a  mechanism  are  described  in  terms  of  the  configuration 
space  of  the  mechanism.  A  proper  decomposition  of  the  configuration  space,  called  a  region  di- 
agram, reflects  the  possible  motions  of  each  object  in  the  mechanism  and  the  relationships 
between  these  motions.  We  show  how  to  obtain  the  region  diagram  of  a  fixed  axis  mechanism  by 
composing  pairwise  configuration  spaces.  The  resulting  region  diagram  can  be  used,  given  se- 
quence of  motions  applied  to  different  objects  in  the  mechanism  and  an  initial  state  of  the 
mechanism,  to  produce  a  state  diagram  describing  the  qualitative  behavior  of  the  mechanism. 


1  This  work  was  supported  by  the  Defense  Advanced  Research  Projects  Agency  under  contract  N00014-85-K- 
0163  from  the  Office  of  Naval  Research  and  by  a  National  Science  Foundation  grant  under  contract  DCR-860375S. 


There  has  been  considerable  attention  in  recent  years  on  representing  and  reasoning  about 
physical  objects  and  processes.  A  problem  that  is  of  relevance  in  this  context  is  explaining  how  a 
given  physical  system  such  as  a  boiler,  an  electrical  circuit  or  a  clutch  work.  Previous  work  [De- 
Kleer and  Brown  85],  [Forbus  85],  [Kuipers  85]  has  focused  on  the  problem  of  deriving  the  func- 
tion of  a  physical  system  from  the  description  of  its  individual  components  and  their  interconnec- 
tions. The  description  of  the  function  of  the  system  is  stated  in  terms  of  possible  qualitatively  dif- 
ferent behaviors  and  their  relationship.  Our  interest  is  in  reasoning  about  assemblies  of  rigid  ob- 
jects, in  particular  in  the  domain  of  mechanical  devices.  A  mechanism  is  a  collection  of  objects 
(also  called  parts)  assembled  so  that  mechanical  forces  can  be  compelled  to  do  work  accompanied 
by  certain  determinate  motions  (Reuleaux's  definition  of  a  machine;  [Reuleaux  1876]).  We  res- 
trict ourselves  to  the  kinematic  analysis  of  mechanisms,  i.e  the  positions  and  possible  motions  of 
parts  in  the  mechanism  and  ignore  dynamics  issues  such  as  frictional  forces,  torques,  etc.  The 
behavior  of  a  mechanism  can  thus  be  described  by  the  possible  motions  that  every  object  can 
have,  and  their  relationships. 

Deducing  the  behavior  of  a  device  from  its  structural  description  is  necessary  to  carry  a 
number  of  tasks  such  as  explanation,  diagnosis  and  design.  In  all  three  cases  we  must  be  able  to 
determine  the  possible  behaviors  of  a  device  so  that  we  can  either  produce  a  causal  explanation, 
determine  the  fault  by  comparing  the  observed  behavior  with  the  desired  behavior  or  introduce 
structural  changes  to  achieve  the  desired  behavior.  For  these  tasks,  envisioning  and  prediction  of 
behavior  are  of  particular  importance.  Given  a  structural  description  of  a  mechanism  in  terms  of 
its  parts,  their  shapes  and  an  initial  position,  envisioning  consists  of  determining  all  the  possible 
motions  of  each  part  and  their  relationships.  Given  a  sequence  of  motions  applied  to  various  parts 
of  the  mechanism,  prediction  of  behavior  consists  of  determining  the  actual  motions  of  each  part 
and  their  relationships  Furthermore,  a  qualitative  description  of  these  behaviors  is  desirable  to 
produce  good  explanations  and  guide  the  diagnosis  or  design  process. 

Previous  work  in  Qualitative  Physics  has  addressed  the  problem  of  determining  the  qualita- 
tive behavior  of  a  physical  system  from  its  structure.  In  De  Kleer  and  Brown's  formalism  [De- 
Kleer  and  Brown  85],  a  device  consists  of  three  types  of  constituents:  materials,  components  and 
conduits.  Components  are  elementary  parts  that  operate  on  and  chr.nge  materials.  Conduits  are 
components  that  do  not  change  materials:  they  transport  the  material  from  one  component  to  the 
other.  Behavior  is  achieved  by  transporting  materials  from  one  component  to  the  other  through 
conduits.  Applying  this  paradigm  to  the  domain  of  mechanical  devices  amounts  to  considering  mo- 
tions as  materials  that  are  transported  by  mechanical  parts  and  that  are  modified  by  pairs  of  parts. 
Note  that  individual  parts  cannot  be  considered  as  primitive  components  since  the  function  of  a 
part  depends  on  the  type  of  contacts  it  has  with  other  parts.  For  example,  a  cylinder  can  be  both 
used  as  a  supporting  axis  for  a  gear  and  as  a  centering  pin.  This  is  in  contrast  with  the  modeling 
of  electrical  devices  as  described  in  [DeKleer  and  Brown  85],  where  electrical  parts  (resistances, 
light  bulbs,  batteries)  are  primitive  components  and  their  function  does  not  change  in  different 
devices  (radios,  heaters,  etc.).  Thus,  for  mechanical  devices  each  type  of  pair  is  a  component  and 
its  kinematic  function  is  described  in  terms  of  the  motions  of  each  part.  For  example,  the  kinemat- 
ic function  of  a  pair  of  gears  is  to  provide  a  constant  ratio  between  the  rotation  of  the  two  gears. 
This  component  approach  has  a  number  of  disadvantages.  First,  it  requires  that  all  the  relations 
between  pairs  of  objects  to  be  defined  in  advance.  But  since  the  number  of  possible  objects  is  very 
large  and  small  changes  in  their  geometry  can  lead  to  radically  different  behaviors,  an  exhaustive 
list  of  such  relations  is  impractical.  Second,  DeKleer  and  Brown's  model  does  not  account  for 
changes  in  the  topology  of  the  device.  Changes  in  topology  appear  very  often  in  mechanical  dev- 
ices and  thus  cannot  be  ignored.  Forbus'  Qualitative  Process  Theory  [Forbus  85]  suffers  from  the 
same  limitations. 

Quantitative  geometrical  reasoning  about  moving  objects  has  been  studied  in  relation  to  the 
motion  planning  problem  [Lozano-Perez  and  Wesley  79],  [Lozano-Perez  82]  and  [Schwartz  and 
Sharir  82].  Given  an  initial  and  a  final  configuration  of  the  objects,  the  goal  is  to  find  a  sequence 
of  motions  in  space  that  describe  the  path  (if  such  a  path  exists)  from  one  configuration  to  the 


other  A  general  but  inefficient  solution  has  been  provided  in  [Schwartz  and  Sharir  82]  and  many 
special  cases,  mainly  of  one  moving  object  amidst  fixed  obstacles,  have  been  analyzed  [Schwartz 
et.  al  87]  Csee  their  bibliography).  These  techniques  are  relevant  since  finding  the  possible 
behaviors  of  the  parts  in  a  mechanism  can  be  viewed  as  finding  the  set  of  all  the  feasible  confi- 
gurations that  a  mechanism  can  have,  and  the  paths  and  motion  sequences  required  to  achieve 
them.  This  approach  however  is  prohibitively  expensive,  even  for  mechanisms  with  a  small 
number  of  moving  parts.  Features  such  as  restricted  number  of  degrees  of  freedom  of  objects  and 
the  small  number  of  possible  contacts  between  objects  are  frequently  found  in  mechanisms  and 
simplify  their  analysis.  Furthermore,  there  is  still  the  need  to  analyze  the  possible  motions 
described  in  quantitative  terms  and  produce  a  qualitative  description  of  them. 

More  recently,  [Fairings  87],  [Nielsen  87],  [Gelsey  87]  and  [Joskowicz  87]  address  the  partic- 
ular problem  of  the  analysis  of  the  kinematic  behavior  of  mechanisms.  [Faltings  87]  describes  a 
method  to  compute  a  qualitative  "place  vocabulary"  that  describes  the  set  of  allowable  object 
placements  in  the  plane.  [Nielsen  87]  provides  a  method  of  using  this  place  vocabulary  to  produce 
a  qualitative  description  of  the  object's  kinematic  behavior.  [Gelsey  87]  describes  an  analysis 
method  based  on  the  composition  of  behavior  of  kinematic  pairs.  [Joskowicz  87]  describes  a 
two-step  method  for  the  qualitative  analysis  of  mechanical  devices.  The  first  step  takes  the 
geometrical  description  of  the  parts  in  their  initial  configuration  and  analyzes  how  two  objects  res- 
trict each  others'  motions.  The  result  is  a  description  of  the  possible  relative  motions  of  the  pair. 
The  second  step  takes  the  pairwise  relative  motions  and  combines  them  to  produce  the  possible 
morions  of  each  part  and  their  relationships  using  a  constraint  propagation  technique.  The  final 
output  is  a  state  diagram  similar  to  the  one  in  [DeKleer  and  Brown  85]  that  describes  the  behavior 
of  the  mechanism.  This  method  is  restricted  to  fixed  axis  mechanisms. 

The  purpose  of  this  report  is  to  establish  a  formal  framework  for  the  kinematic  analysis  of 
mechanisms  based  on  the  idea  of  kinematic  pairs  and  their  configuration  spaces.  Using  the  confi- 
guration space  approach,  we  provide  a  justification  for  the  two-step  method  in  [Joskowicz  87]  and 
show  how  it  is  po^ible  to  obtain  the  configuration  space  of  the  mechanism  from  the  composition 
of  local  configuration  spaces  of  pairs  of  objects.  We  also  show  how  this  composition  can  be  effec- 
tively carried  for  fixed  axis  mechanisms.  We  do  not  deal  with  the  issue  of  how  actually  to  obtain 
the  local  configuration  spaces  from  the  geometry  of  the  objects.  We  assume  throughout  the  report 
that  these  configuration  spaces  are  given. 

This  report  is  organized  as  follows:  the  first  section  introduces  the  concept  of  configuration 
spaces  in  the  context  of  mechanical  devices  and  shows  how  it  is  possible  to  obtain  the  configura- 
tion space  of  the  overall  mechanism  by  composing  pairwise  configuration  spaces.  The  relationship 
between  subsets  of  the  configuration  space  and  the  possible  motions  of  parts  is  explained  and  a  set 
of  possible  motions  is  introduced.  The  second  section  defines  a  decomposition  of  the  configura- 
tion space  of  the  mechanism,  called  the  region  diagram.  We  then  show  how  to  obtain  the  region 
diagram  for  the  whole  mechanism  by  composing  region  diagrams  of  pairs  of  objects  and  give  an 
algorithm  to  do  this  composition  for  fixed  axis  mechanisms.  The  third  section  presents  an  algo- 
rithm to  construct  the  state  diagram  from  a  region  diagram,  a  sequence  of  input  motions  and  an 
initial  state  of  the  mechanism.  The  states  in  the  diagram  correspond  to  qualitatively  different  mo- 
tions of  the  parts.  We  conclude  by  examining  the  difficulties  of  extending  this  approach  to  mov- 
able axis  mechanisms. 


1.    Configuration  Spaces 

The  concept  of  configuration  spaces  was  initially  introduced  in  the  context  of  motion  plan- 
ning [Lozano-Perez  and  Wesley  79],  [Lozano-Perez  82],  [Brooks  and  Lozano-Perez  83].  Planning 
the  motion  of  an  object  amidst  obstacles  requires  the  computation  of  the  free  placements  of  the 
object  to  determine  whether  it  is  possible  to  move  from  one  location  to  another  and  what  path 
must  the  object  follow  in  order  to  achieve  the  desired  final  position.  Since  a  three-dimensional 
object  can  have  6  degrees  of  freedom  in  space  (3  rotations  and  3  translations),  its  placement  (posi- 
tion and  otientation)  in  space  can  be  uniquely  characterized  with  respect  to  its  reference  point  by 
a  point  in  six  dimensional  space.  [Lozano-Perez  and  Wesley  79]  show  that  the  planning  the  motion 
of  a  rigid  object  amidst  obstacles  amounts  to  finding  a  path  for  the  reference  point  of  the  moving 
object  in  a  6-dimensional  space  amidst  6-dimensional  obstacles.  Planning  the  motion  of  n  objects 
can  be  viewed  as  planning  the  motion  of  a  point  in  6n-dimensional  space  amidst  6n  dimensional 
obstacles.  The  configuration  space  is  the  set  of  free  placements  of  the  point  representing  the 
objects  in  the  new  space.  It  is  this  set  that  determines  which  motions  are  feasible. 

Configuration  spaces  provide  a  powerful  abstraction  mechanism  since  assemblies  of  objects 
having  different  geometries  can  have  the  same  configuration  space,  and  the  nature  of  this  confi- 
guration space  has  a  direct  correspondence  to  their  possible  motions.  Consider  the  three  mechan- 
isms in  Figure  1.1  and  assume  that  the  base  is  fixed.  The  configuration  space  of  the  slider  with 
respect  to  the  base  is  the  same  in  the  three  cases.  It  indicates  that  the  only  possible  free  place- 
ments of  the  slider  with  respect  to  the  base  are  along  a  line.  Since  all  the  possible  placements  of 
the  slider's  reference  point  lie  on  a  line,  the  only  possible  motion  of  the  slider  is  a  translation 
along  a  fixed  axis  parallel  to  O. 

Different  topologies  of  the  configuration  space  correspond  to  different  types  of  possible 
motions.  In  the  previous  example  the  hyperplane  defined  by  the  equations  in  Figure  1.1 
corresponds  to  rectilinear  paths  and  thus  to  translational  motions.  The  computation  of  the  confi- 
guration space  provides  us  with  a  mapping  from  the  geometrical  description  of  the  objects  to  the 
set  of  positions  of  their  reference  points  in  a  manifold.  The  nature  of  this  manifold,  in  turn, 
corresponds  to  the  types  of  motions  that  objects  can  have.  This  mapping  provides  us  with  an 
abstraction  mechanism  that  starts  with  a  geometrical  description  of  a  set  of  objects  and  produces  a 
"functional"  description.  As  argued  in  [Hopcroft  and  Krafft  87],  such  a  description  is  desirable 
since  it  can  support  symbolic  deductions  not  based  on  their  actual  geometry.  The  mechanisms  in 
Figure  1.1  all  belong  to  the  prismatic  joint  class  and  we  can  associate  with  this  class  of  pairs  a  set 
of  symbolic  properties  that  can  be  used  in  diagnosis  or  design  tasks. 

Although  it  is  possible  in  principle  to  compute  the  configuration  space  of  the  overall  set  of 
objects  at  once  as  it  is  done  in  the  motion  planning  problem,  this  approach  is  not  practical  for 
mechanisms.  First,  the  dimensionality  of  the  configuration  space  grows  linearly  with  the  number 
of  parts  in  the  mechanism  since  there  are  6  potential  degrees  of  freedom  for  each  object. 
Nonetheless,  most  mechanisms  have  very  few  degrees  of  freedom  (one  or  two)  and  thus  it  is  pos- 
sible to  analyze  a  projection  of  their  configuration  space  into  a  lower  dimensional  space.  Second, 
the  concept  of  locality  of  function  is  not  exploited  by  this  approach.  In  many  mechanisms, 
subassemblies  are  designed  to  perform  a  specific  kinematic  task  and  can  therefore  be  analyzed 
independently.  Also,  the  number  of  possible  interactions  between  parts  is  small  since  every  part  is 
not  in  contact  with  every  other  part.  Finally,  most  parts  are  designed  to  have  only  simple  possible 
motions  such  as  rotation  or  translation. 

These  considerations  lead  us  to  consider  an  approach  in  which  relations  between  pairs  of 
objects  are  examined  first  and  then  put  together  to  obtain  the  overall  configuration  space.  The 
advantage  is  that  we  have  to  deal  only  with  pairwise  configuration  spaces  of  fixed  dimensionality 
(12  dimensional).  We  will  show  in  this  section  that  the  overall  configuration  space  of  the  mechan- 
ism can  be  computed  exactly  by  composing  pairwise  configuration  spaces. 


1.1.  Objects 

For  the  purposes  of  this  formal  framework,  objects  will  be  considered  as  infinite  sets  of 
points  in  three  dimensional  space.  The  only  requirement  we  have  is  that  each  set  be  connected  and 
describable  by  a  finite  set  of  algebraic  relations.  The  placement  of  an  object  in  space  is  the 
description  of  its  position  and  orientation. 

Definition  1.1.1:  A  three  dimensional  rigid  object  is  a  connected,  non-empty  point  set 
A  C  R\ 

Definition  1.1.2:  The  function  "interior(A)"  maps  a  point  set  A  onto  its  interior,  such  that  a 
point  p  is  in  interior(A)  if  the  neighborhood  of  p  in  A  is  homeomorphic  to  Rk  (0  s  k  s  3). 
Therefore  intenor(A)  C  A. 

Definition  1.1.3:  The  function  'boundary(A)"  maps  a  point  set  onto  its  boundary  such  that 
boundary(A)  =  A  -  intenor(A).  In  a  mechanism,  we  require  that  no  two  objects  overlap,  i.e. 
that  their  interiors  are  always  disjoint,  (intenor(A)  n  interior(B)  =  0).  The  predicate 
"overlap(A,  B)"  is  true  iff  A  and  B  are  not  disjoint. 

Definition  1.1.4:  Two  objects  are  in  contact  iff  their  boundaries  intersect: 
contact(A,  B)  =>  A  n  B  =  boundary(A)  n  boundary(B)  *  0. 

Definition  1.1.5:  A  mechanism  M  is  a  collection  of  n  three  dimensional  objects,  M  =  {Ai} 
such  that  at  least  one  object  is  fixed  in  space.  This  fixed  object  is  called  the  supporting  frame  of 
the  mechanism. 

1.2.  Configuration  Space 

Let  M  be  a  mechanism  consisting  of  n  three-dimensional  objects  Al,  ...  An.  Let  F  be  a 
reference  frame  of  the  Euclidean  space  R3  with  origin  at  O.  Let  each  object  Ai  have  a  reference 
point  O     with  respect  to  O. 

In  three  dimensional  space,  each  object  has  6  degrees  of  freedom,  and  its  position  and  orien- 
tation can  be  uniquely  characterized  with  respect  to  O^  in  R3  by  a  6-tuple  <x|>,  1  £  j  ^  6 
where  x  ,  x,,  x3  represent  the  position  of  the  reference  point  in  Cartesian  coordinates,  and  xj4, 
x;,  x6  (also  denoted  8  ,  <J>  and  4»,)  specify  the  orientation  of  Ai  with  respect  to  the  reference 
frame  F  (see  Figure  1.2.) 

Definition  1.2.1:  Cspace(Ai)  is  the  set  of  all  possible  placements  in  six  dimensional  space  of 
Ai  and  is  defined  by  the  parameters  x  ,  ....  x  .  Cspace(M)  is  the  6n  dimensional  space 
corresponding  to  a  mechanism  with  n  objects. 

Definition  1.2.2:  Let  p  be  a  point  in  Cspace(Ai).  Let  placement(Ai,  p)  be  a  family  of  func- 
tions 

placement(Ai,  p)  :  Cspace(Ai)  =>  2Ri 

where  2R  is  the  set  of  subsets  of  R\  such  that  for  a  placement  p  =  <*,,,  ...,  X|6>  €  Cspace(Ai), 
placement(Ai,  p)  is  the  set  of  points  occupied  by  Ai  when  placed  such  that  O  is  at  position  <*,,, 
x,,  x,>  and  its  orientation  is  given  by  <x  .,  x  .,  x  ,> .  This  definition  is  also  extended  to  sets  of 
objects,  i.e  placement(M,  p)  represents  the  set  of  points  occupied  by  the  objects  in  M  when  placed 


m  space. 

Definition  1.2.3:  Let  p  be  a  point  in  Cspace(M).  Let  "projection(Ai,  p)"  be  the  family  of 
functions 

projection(Ai,  p):  Cspace(M)  =>  2c^iceiA,) 
such  that  if  p  €  Cspace(M)  and  q  €  Cspace(Ai),  q  =  projection(Ai,  p)  iff  p  =   <*u,  ••••,  *n,  •••. 
x   ,  ...,  x  ,>  and  q  =   <x  .,  ...,  x6>.  The  function  "projection(Ai,  p)'-  takes  a  point  in  Cspace(M) 
and  projects  it  into  Cspace(Ai). 

Let  p  be  a  point  in  Cspace(M)  and  pi  its  projection  in  Cspace(Ai),  p  =  projection(Ai,  p).  To 
simplify  the  notation,  we  will  write  placement(Ai,  p)  to  mean  placement(Ai,  projection(Ai,  p)). 

Definition  1.2.4:  The  configuration  space  of  M  with  respect  to  a  fixed  reference  frame  F, 
CO(M)  is  a  subset  of  Cspace(M)  such  that  for  every  placement  of  M.  no  two  objects  Ai,  Aj  over- 
lap: 

CO(M)  =  {p  €  Cspace(M)  j  not(overlap(placement(Ai,  p),  placement(Aj,  p)),  i  *■  j}, 

The  set  CO(M)  represents  the  set  of  a     free  placements  of  the  parts  of  M. 

Definition  1.2.5:  The  configuration  space  COM(Ai)  of  an  individual  object  with  respect  to  M 
is  the  set  of  possible  placements  of  Ai. 

COM(Ai)  =  {projection(Ai,  p)  \  p  €  CO(M)} 

Note  that  since  we  require  the  mechanism  to  have  a  fixed  supporting  frame,  (say  Ak),  the  place- 
ment of  Ak  can  be  characterized  by  a  single  point  in  its  configuration  space,  i.e  COM(Ak)  =  {p} 
where  p  is  a  point  of  Cspace(Ak).  The  requirement  of  a  fixed  supporting  frame  is  necessary,  for 
otherwise  COM(Ai)  will  always  be  equal  to  Cspace(Ai)  since  it  is  always  possible  to  move  all  the 
parts  of  a  mechanism  together  so  that  a  part  will  achieve  any  position  and  orientation  in  space. 
This  is  no  longer  possible  if  at  least  one  part  is  fixed  in  space  and  creates  an  obstacle. 

Notation:  Note  that  CO  (M)  =  CO(M).  For  convenience  we  will  denote  the  configuration  space 
of  a  pair  of  objects  as  CO(A,  B)  instead  of  CO({A,  B}). 

The  set  CO(M)  is  a  6n  dimensional  manifold  that  can  be  characterized  by  a  finite  set  of  alge- 
braic relations,  provided  that  the  objects  are  also  described  by  a  set  of  algebraic  relations.  See 
[Schwartz  and  Sharir  82]  for  a  complete  algebraic  characterization  of  the  configuration  space. 

1.3.    The  composition  operator 

We  now  show  how  to  compute  the  configuration  space  of  the  overall  mechanism  by  compos- 
ing configuration  spaces  of  subassemblies.  In  the  following  we  assume  that  the  mechanism  M  con- 
sists of  at  least  two  objects. 

Definition  1.3.1:  Let  M  be  a  mechanism  and  Ml  and  M2  two  nonempty  submechanisms.  Let 
Cspace(M),  Cspace(Ml)  and  Cspace(M2)  be  their  corresponding  parameter  spaces.  Let 
Cspace(M)  be  defined  by  the  parameters  {x ..},  lsis  |M|,  Is  j  s  6,  Cspace(Ml)  be  defined  by 
{yk.},  1  s  k  s  |M1|  and  Cspace(M2)  be  defined  by  {z,!,  1  s  1  s  |M2|.  Clearly,  {ykj}  C  {xj  and 
{zjcfx,}. 

We  define  the  composition  operation,  denoted  |x  j  as  follows: 


Definition  1.3.2:  The  parameter  space  of  the  composition  of  Ml  and  M2,  denoted  Cspace(Ml 
|x  |  M2)  is  a  6m  dimensional  space,  where  m  =  .Ml  U  M2 ,  and  is  defined  by  the  parameters  {yk } 
U  {Z]j}. 

Definition  1.3.3:  Let  SI  C  Cspace(Ml)  and  S2  C  Cspace(M2).  Then 

SI  |x  |  S2  =  {p  €  Cspace(Ml  |x  |  M2)  |  projection(Ml,  p)  (.  SI  a  projection(M2,  p)  (.  S2} 

Note  that  the  composition  operation  is  very  similar  to  the  natural  join  operation  used  in  rela- 
tional calculus,  although  the  composition  operation  is  defined  over  infinite  sets. 

1.4.    Properties  of  the  composition  operation 

Here  are  some  useful  properties  of  the  composition  operation.  Let  S(M1),  S'(M1),  S(M2) 
and  S'(M2)  be  distinct  subsets  of  Cspace(Ml)  and  Cspace(M2)  respectively. 

1.  If  Ml  and  M2  have  no  parts  in  common  (i.e.  they  are  disjoint,  the  composition  of  two  sub- 
sets of  Cspace(Ml)  and  Cspace(M2)  is  their  cross  product: 

Ml    n    M2  =     0      =      S(M1)  |x|   S(M2)    =    S(M1)  x  S(M2) 

2.  The  composition  of  a  subset  of  Cspace(M)  with  itself  is  the  same  subset: 
S(M)  |x|S(M)  =  S(M) 

3.  The  composition  of  two  subsets  of  Cspace(M)  is  equal  to  their  intersection: 
S(M)  !x|S'(M)  =  S(M)  n  S'(M) 

4.  The  composition  operation  is  symmetric: 
S(M1)  |x  |  S(M2)  =  S(M2)  |x  |  S(M1) 

5.  The  composition  operation  is  associative: 

S(M1)  |x  |  (S(M2)  |x  |  S(M3))  =  (S(M1)  |x  I  S(M2))  |x  |  S(M3) 

6.  The  composition  operation  is  distributive  with  respect  to  intersection  and  union: 
(S(M1)  U  S'(M1))  |x  |  S(M2)  =  (S(M1)  |x  |  S(M2))  U  (S'(M1)  |x  |  S(M2)) 
(S(M1)  D  S'(M1))  |x|S(M2)  =  (S(M1)  jx|S(M2))  n  (S'(M1)  |x  j  S(M2)) 

7.  The  composition  of  a  subset  of  Cspace  and  Cspace  is  the  subset: 
Cspace(M)  |x|S(M)  =  S(M) 

8.  The  composition  operation  is  monotonic: 

S(M1)  C  S'(M1)    =>  S(M1)  |x  |  S(M2)  C  S'(M1)  |x  |  S(M2) 

1.4.1.    Computing  CO(M)  from  subassemblies 

Let  M  be  a  mechanism  with  n  parts  (n  >  1)  and  CO(M)  its  configuration  space.  Let  Ml,  ... 
Mn  be  all  the  subassemblies  of  n  -  1  parts  of  M  and  CO(Ml),  ...  CO(.Mn)  their  corresponding 
configuration  spaces  independent  of  M  (i.e  not  COM(Mi)). 

Theorem  1:  The  configuration  space  of  M  can  be  computed  by  the  composition  of  indepen- 
dent configuration  spaces  of  all  n(n  -  l)/2  pairs  of  objects  of  M: 

CO(M)    =     CO(Al,  A2)    |x|   CO(Al,A3)    |x  |   ...    |x  |   CO(An-l,  An) 

Theorem  2:  The  configuration  space  of  M  can  be  computed  by  the  composition  of  indepen- 
dent configuration  spaces  of  all  n  subassemblies  of  M  containing  n  -  1  objects: 

CO(M)      =      CO(Ml)    |x|  CO(M2)    |x|   ...    |x  |  CO(.Mn) 
Proof  of  Theorem  1:  By  the  definition  of  configuration  space1: 


For  compactness  \*e  write  (Ai)    to  mean  placement Ai.  p).  the  position  of  Ai  in  Euclidean  space  for  a  given 


CO(M)  =  {p  €  Cspace(M)  I  not(overlap((Ai)_,  (Aj)  )  V  i,j    1  *  \\ 

P  P 

=  {p  k  Cspace(M)  i        not(overlap((Al)   ,     (A2)  )        and 

p  p 

not(overlap((Al)   ,     (A3)   )        and 

p  p 

not(overlap((An-l)p,  (An)p)} 

Oq  the  other  hand, 

CO(Al,  A2)      =  {p  €  Cspace(Al,A2)  |  not(overlap((Al)  ,  (A2)  )} 
CO(Al,  A3)      =  {p  €  Cspace(Al,A3)  |  not(overlap((Al)  ,  (A3)  )} 

CO(An-l,  An)  =  {p  €  Cspace(An-l,  An)  |  not(overlap((An-l)  ,  (An)  )} 

Let  xi  be  the  parameter  vector  corresponding  to  object  Ai,  xi  =  <x  ,  ...  x.6>.  Bv  the  definition 
of  the  composition  operation: 

CO(Al,A2)    |x|   CO(Al.  A3)    |x|   ...    |x  |     CO(An-l,  An)  = 
{p  =  <xl,  ...,  xn>   €  Cspace(M)  | 

not(overlap((Al),u  x2>,  (A2)<u  0>)  and 

not(overlap((Al)<Klo:>,  (A3)<u  0>)  and 

not(overlap((An-lK     .      .  ,  (An),.     ,      .J} 

Note  that  by  substituting  <xi,  xj>  by  <xl,  ...,  xn>  the  set  remains  unchanged.  □ 

Proof  of  Theorem  2:  We  will  prove  that  this  theorem  is  equivalent  to  Theorem  1. 

Let  (Ai,  Aj)  be  all  the  n(n  -    l)/2  possible  pairs  of  distinct  parts  of  M  and  let  CO(Ai,  Aj) 
their  corresponding  independent  configuration  spaces. 

We  first  prove  that  Theorem  2  implies  Theorem  1.  The  proof  is  by  recursion  on  the  number  of 
parts  of  M.  Assume  Theorem  2  is  true  for  a  mechanism  of  n  objects.  Then  it  is  also  true  for 
mechanisms  with  n  -  1,  n  -  2,  ...  2  objects.  Apply  the  formula  (2) 

CO(M)  =  CO(Ml)    |x|  CO(M2)    |x|   ...    |x|   CO(Mn) 

to  Ml,  ....  Mn  recursively  n  -  2  times  until  the  number  of  parts  in  each  subassembly  Mi  is  2.  We 
get  n!  terms  of  the  form  CO(Ai,  Aj)  where  1  £  i,j  s  n.  Since  there  are  only  n(n  -  l)/2  different 
pairs,  each  term  is  repeated  2(n  -  2)!  times.  By  using  properties  4  and  5,  we  can  group  together 
all  the  terms  that  are  equal  and  simplify  all  the  compositions  of  the  form  CO(Ai,  Aj)  |x|  CO(Ai, 
Aj)  to  CO(Ai,  Aj)  using  Property  2.  We  are  left  with  n(n  -  l)/2  terms  of  the  form  CO(Ai,  Aj) 
for  all  i,  j  such  that  1  s  i  <  j  s  n: 

CO(M)  =  CO(Al,  A2)    |x|  CO(Al,  A3)    |x  |   ...    |x  |   CO(An-l,  An) 

D 

Now  we  prove  that  Theorem  1  implies  Theorem  2.  We  prove  it  by  induction  on  the  number  of 
parts  of  M:  for  a  mechanism  of  3  parts,  equation  (2)  becomes 

CO(M)  =  CO(Al,  A2)    |x|   CO(Al,  A3)    |x  |  CO(A2,  A3) 

since  there  are  only  three  combinations  of  size  2  of  M  =  {AI,  A2,  A3}.  This  is  Identical  to  (1) 
for  n  =  3.  Now  assume  that  Theorem  2  is  true  for  mechanisms  of  n  parts,  i.e 

CO(M)  =  CO(Ml)    |x|   CO(M2)    |x|   ...    |xjCO(Mn) 

CO(Al,A2)    |x|   CO(Al,  A3)    |x|...    |x|   CO(An-l,  An) 


point  p  in  CO(M) 


Let  M  be  a  mechanism  of  a-  1  parts.  By  Theorem  1, 

CO(M)  =  CO(Al,  A2)    |x|CO(Al,  A3)    |x  |   ...    |x  |   CO(An,  An-  I) 

By  repeating  n!  times  each  CO(Ai,  Aj)  using  property  2,  and  grouping  the  pairwise  configuration 
spaces  into  groups  of  n  parts  using  properties  4  and  5  we  get: 

CO(M)  =  [CO(A2,  A3)    |x|    ...    |x  |   CO(An,  An-1)]  (exclude  Al) 

|x| 

[CO(Al,  A3)    |x|    ...    |x|   CO(An,  An-1)]  (exclude  A2) 


[CO(Al,  A2)    |x|    ...    |x|   CO(An-l,  An)]  (exclude  An- 1) 

By  the  induction  hypothesis,  for  mechanisms  Mi  =  M  -  {Ai}  of  n  objects,  1  s  i  a  n  +  1: 

CO(Ml)  =         CO(A2,  A3)    |x|    ...    |x|   CO(An,  An-1) 

CO(M2)  =         CO(Al,  A3)    |x  |    ...    |x  |   CO(An,  An-1) 

CO(Mn-l)    =CO(Al,A2)    |x  |    ...    |x|   CO(An-l,  An) 

substituting  each  CO(Mi)  in  the  previous  expression  leads  to  Theorem  2.  Therefore,  the  confi- 
guration space  of  m  can  be  computed  by  composing  pairwise  configuration  spaces. 

We  now  turn  our  attention  to  mechanisms  in  which  certain  objects  never  come  into  contact. 

Theorem  3:  Let  M  be  a  mechanism  consisting  of  three  objects,  A,  B  and  C  such  that  A  and  C  can 
never  come  into  contact,  i.e  V  p  €  CO(M),  (A)  n  (C)  =  0.  The  configuration  space  of  M  is 
the  composition  of  the  independent  configuration  spaces  of  the  pairs  A,  B  and  B,  C: 

V  p  €  CO(M),  (A)    n  (C)p  =  0     =>     CO(M)  =  CO(A.  B)     |x  |    CO(B,  C) 

Proof:  Since  V  p  €  CO(M),  (A)    n  (C)    =0  implies  that 
V  p  €  CO(M),  not(overlap((A)p  (~l  (C)  ). 
By  the  definition  of  CO(M), 

CO(M)  =  {p  €  Cspace(M)  I       not(overlap((A)p,  (B)p)    and 

not(overlap((B)p,  (C)p)    and 
not(overlap((A)Pp,  (C)p)} 

By  the    definition  of  the  composition  of  CO(A,  B)  and  CO(B,  C): 

CO(A,  B)  |x|CO(B,  C)  = 

{p  €  Cspace(M)  |  not(overlap((A)p,  (B)p)  a  not(overlap((B)p,  (C)p)} 

Thus, 

CO(A,  B)  |x  |  CO(B,  C)  and  V  p  €  CO(M),  not(overlap((A)p  n  (C)p)    =  CO(M) 

Theorem  3  can  be  generalized  to  mechanisms  with  several  parts,  for  not(overlap(Ai)p, 
(Aj)  )  is  identically  true  in  CO(M)  if  Ai  and  Aj  never  come  in  contact: 

F 

Theorem  4:  Let  M  be  a  mechanism  and  Ml  ...  Mk  subassemblies,  not  necessarily  containing 
the  same  number  of  parts  such  that  M  =  Ml  U  ...U  Mk  and  such  that,  if  two  parts  ever  come  in 
contact,  then  there  is  a  subassembly  containing  both  of  them.  Then  the  configuration  space  of  M 
is  the  composition  of  the  configuration  spaces  CO(Ml),  ...  ,  CO(Mk). 

If  V  Ai,  Aj  in  M: 


(3  p  i  CO(M)  (Ai)    H  (Aj)J  *  0  =-  (3  Ml  such  that  Ai,  Aj  €  Ml) 

Then,  CO(M)  =  CO(Ml)    jx  |   CO(M2)    |x|   ...    |x|   CO(Mk) 
Proof:  Follows  from  the  use  of  Theorems  1,  2  and  3. 

Note:  The  fact  that  CO(M)  can  be  computed  as  a  composition  of  n(n  -  l)/2  painvise  configura- 
tion spaces  does  not  mean  that  CO(.M)  can  be  computed  in  polynomial  time,  since  CO(M)  can 
require  an  exponential  number  functions  to  be  described  (exponential  in  the  number  of  degrees  of 
freedom  of  the  objects).  The  composition  operation  is  difficult  to  evaluate  in  the  general  case,  but 
we  will  see  that  in  some  particular  situations  it  can  be  performed  efficiently. 

1.5.  Representations  of  the  configuration  space 

The  standard  way  of  characterizing  a  configuration  space  is  by  describing  it  with  a  set  of 
algebraic  inequalities  that  define  the  portion  of  space  corresponding  to  free  placements.  For  exam- 
ple, in  the  case  of  the  prismatic  joint,  the  set  of  free  positions  of  A  is  described  by  the  following 
equations.  Let  p  =  <x,  y,  z,  6,  <J>,  ij>>  be  a  point  of  the  configuration  space  CO    (A).  Then, 

V  p  €  COM(A)    : 

PI:      ax  -  by  -  cz  =  d 

P2:      a'x  -  b'y  -  c'z  =  d' 

0  =  constant,      <j>  =  constant,      <li  =  constant 

x  €  [Xmin,  Xmax],      y  €  [Ymin,  Ymax],      z  €  [Zmin,  Zmax] 

that  is,  the  three  coordinates  that  represent  the  position  of  A  in  space  lie  on  a  line  (the  intersec- 
tion of  two  planes)  and  the  three  coordinates  that  correspond  to  the  orientation  of  A  are  constant. 

We  can  replace  the  algebraic  description  of  a  subset  S  of  COM(Ai)  by  a  parameterized 
description  with  a  different  reference  frame  than  the  original  one.  For  example,  in  the  case  of  a 
prismatic  joint,  we  can  describe  the  configuration  space  by  a  new  set  of  algebraic  inequalities 
using  a  rigid  mapping  to  a  new  coordinate  frame  in  which  the  translation  axis  T  is  parallel  to  the 
new  x'  axis,  so  that  the  only  varying  parameter  is  along  that  axis,  the  axis  of  translation.  We 
define  therefore  a  fixed  rigid  mapping  a  between  the  points  in  the  space  defined  by  <x,  y,  z,  8, 
<(>,  <|»>  and  the  points  in  the  new  space  defined  by  <x',  y',  z',  6',  4>',  i|T>.  In  the  new  coordinate 
frame  for  every  p'  =  <x\  y',  z',  6',  <J>',  <\i'>, 

y'  =  constant,      z'  =  constant,      0'  =  constant,   <J>'  =  constant,   t'  =  constant 
x'  €  [Xmin',  Xmax] 

which  we  describe  by  the  parameterized  set  descriptor  p-translationV((Ai,  O,  x'),  x'  €  [Xmin', 
Xmax'].  The  new  configuration  space  of  Ai  is  COM'(Ai)  =  {cr(p)  |  p  €  COM(Ai)}. 

As  we  will  see  in  the  next  chapter,  this  representation  is  more  convenient  to  do  the  composi- 
tion of  pairwise  configuration  spaces. 

1.6.  Motions 

As  noted  in  the  introductory  example  (Figure  1.1),  there  is  a  direct  relationship  between  the 
nature  of  the  configuration  space  and  the  possible  motions  of  objects.  In  this  section,  we  define 
this  relationship  and  give  a  characterization  of  the  spaces  in  which  simple  motions,  such  as  rota- 
tion or  translation  can  occur. 

Definition  1.6.1:  A  motion  in  CO(M)  is  a  real  valued,  continuous  function  ir(M,  p  ,  p.):  [0, 
1]  -  CO(M)  that  defines  the  trajectory  from  an  initial  placement  pQ  to  a  final  placement  p,  in 
CO(M). 


10 


Definition  1.6.2:  Let  path(-(M.  p.,  p..))  be  the  image  of  the  function  tt(M,  p  ,  p.),  that  is, 
the  subset  of  all  points  of  CO(M)  that  are  traversed  by  the  motion  of  M  from  the  initial  to  the 
final  placement.    We  call  this  set  the  path  of  M  with  respect  to  the  motion  tt(M.  p0,  p.). 

Definition  1.6.3:  The  predicate  "exists-path(Vf ,  pQ,  p.)"  is  true  iff  there  exist  a  motion  -rr  of 
M  from  p0  to  p.: 

exists-path(M,  pQ,  p..)  =  true        iff  3  it  such  that  tt(M,  p  ,  p..)  is  a  motion 

false       else 

Given  a  path  of  M,  path(-(M,  p0,  p.)),  we  can  define  the  path  of  an  individual  part  Ai  as  being 
the  set  of  points  in  COM(Ai)  traversed  by  Ai, 

path(-rr(Ai.  p0,  p,))  =  {projection(Ai,  p)  I  p  €  path(-rr(M,  py  p.))  } 

Definition  1.6.4:  A  path  tt  of  Ai  is  a    rectilinear  path  if  all  the  points  p   =    <x,  y,  z,  6,  (J> , 

ij/>  in  path(-rr(Ai,  p0,  p..))  satisfy  the  following: 

PI:       ax  ■+■  by  -  cz  =  d 

P2:       a'x  -  b'y  -  c'z  =   d' 

9  =  constant,      <i>  =  cons;ant,      t  =  constant 

x  €  [Xmin,  Xmax],      y  €  [Ymin,  Ymax],      z  Z  [Zmin,  Zmax] 

that  is,  the  three  coordinates  of  the  points  in  the  path  that  correspond  to  the  position  of  the  object 
lie  on  a  line  (the  intersection  of  two  non-parallel  planes  PI  and  P2)  and  three  coordinates  that 
correspond  to  the  orientation  of  the  object  are  constant. 

Definition  1.6.5:  The  path  of  Ai  is  a  circular  path  if  all  the  points  p  in  path(iT(Ai,  p0,  pf))  lie 
on  the  surface  described  below.  Let  p  =  <x,  y,  z,  9,  <J>,  \l>>  be  a  point  in  the  path.  Then  the  fol- 
lowing relations  are  satisfied: 

SI:       (x  -  x0):  -  (y  -  yQ)2  -  (z  -  zQ):  =  R2 
PI:       ax  *  by  +  cz  =  d 
(appropriate  relations  for  9,  <t>  and  t) 

that  is,  all  the  three  coordinates  of  the  points  in  the  path  that  correspond  to  the  position  of  the 
object  lie  on  a  circle  (the  intersection  of  a  sphere  SI  of  radius  R  and  of  center  (x0,  y0,  zQ)  and  a 
plane  PI). 

We  now  turn  our  attention  to  connected  subsets  of  the  configuration  space  of  Ai  where  only 
specific  paths  can  appear. 

Definition  1.6.6:  A  possible  path  of  Ai  is  a  subset  of  COV((Ai). 

Definition  1.6.7:  Let  CO  (Ai)  be  the  configuration  space  of  Ai  and  let  S  be  a  connected  sub- 
set of  it.  Let  p  and  p'  be  two  distinct  points  in  S.  We  say  that  the  only  possible  path  in  S  is  a 
translation  iff  for  all  possible  paths  ir(Ai,  p,  p')  in  S,  path(ir(Ai,  p,  p'))  is  a  rectilinear  path. 
Furthermore,  there  exist  two  points  in  S,  p0  and  pf  such  that  S  =  path(ir(Ai,  p0,  p.))  and 
path(-rr(Ai,  pQ,  pf))  is  a  rectilinear  path. 

Let  0A;  be  the  axis  of  the  translation  (defined  by  the  intersection  of  the  two  planes  PI  and 
P2  above).  We  will  rewrite  S  in  this  case  as 

p-translation,,(Ai,  O,  ,  X,  ),  for  X,    €  [Xmin,  Xmax] 
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where  p-translation  ( Ai,  O  ,  X  )  is  a  parametrized  set  descriptor  that  stands  for  rectilinear 
paths  (possible  translations)  along  axis  C»Ai  by  X  ..  This  parameterized  set  descriptor  corresponds 
to  the  mapping 

S  =  CT(path(iT(Ai,  p  ,  p..))  =  p-translationM(Ai,  O    ,  XAj),  for  XAi  €  [Xmin,  Xmax] 

where  a  is  a  fixed  mapping  to  a  new  coordinate  frame.  This  notation  captures  the  fact  that  in  the 
region  S  of  COV((Ai),  the  only  possible  paths  that  Ai  can  follow  are  rectilinear  paths  parallel  to 
axis  O ,  . 

We  can  define  similar  parameterized  set  descriptors  for  paths  defining  possible  motions  such  as 
rotations,  independent  translation  and  rotation  along  the  same  axis  and  dependent  translation  and 
rotation  along  the  same  axis  (see  Figure  1.3).  Following  is  the  list  of  such  descriptors: 

1.  Rotation 

p-rotationM(Ai,  0Aj,  9^),  8Aj  €  [6min.  6max] 
for  possible  rotations  around  axis  O     by  8Aj. 

2.  Translation 

p-translation    (Ai,  0Ai,  XA|),  XAi  €  [Xmin,  Xmax] 
for  possible  translations  along  axis  O     by  XA  . 

3.  Cylindrical  Motion 

p-cylinderM(Ai,  Ov,  XAl,  8Ai),  XAj  €  [Xmin,  Xmax]  and  6A|  €  [9min,  9max] 
for  possible  independent  translation  and  rotation  along  axis  Oaj. 

4.  Helical  Motion 

p-helicalM(Ai,  0A|>  XA|,  8A|),  XAj  =  f(9A|),  9A|  €  [9min,  9max] 

for  possible  helical  motions  (translation  dependent  on  rotation)  along  axis  0A(. 

5.  No  motion 

fixed(Ai,  p) 

for  no  motion  of  Ai  with  respect  to  the  reference  frame.  The  point  p  is  such  that  {p}  is  a 
connected  component  of  COM(Ai). 
We  will  refer  to  this  set  as  the  set  of  possible  simple  motions. 
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2.     Region  Diagrams 

Now  that  we  have  characterized  precisely  the  configuration  space  of  a  mechanism,  we  are 
interested  in  partitioning  it  into  regions  that  will  reflect  qualitatively  different  possible  motions  of 
the  different  parts  of  the  mechanism  and  their  relationships.  We  first  note  that  there  is  a  direct 
relationship  between  subsets  of  the  configuration  space  and  possible  motions.  Consider  the  exam- 
ple in  Figure  2.1.  This  mechanism  consists  of  two  parts,  the  slider  A  and  the  guide  B.  Assuming 
that  the  mechanism  lies  on  the  X-Y  plane,  its  configuration  space  is  best  described  as  a  union  of 
two  subsets,  each  described  by  a  parametrized  set  descriptor  that  correspond  to  the  possible 
motions  of  A  and  B": 

CO(A,  B)    =    Rl    U    R2 

Rl:       p-translation(A,  Ol,  XA)  x  p-translation(B,  Ol,  Xa) 

K  ~  xbI  s  »i 

R2:      p-translation(A,  02,  Y  J  x  p-translation(B,  02,  YB) 

l*A  "  YbI  £  ': 
Note  that  this  decomposition  of  CO(A,  B)  into  two  regions  reflects  the  fact  that  there  are  two 
qualitatively  different  possible  motions  of  A  and  B:  possible  translation  along  axis  Ol  and  possible 
translation  along  axis  02. 

The  purpose  of  this  section  is  to  characterize  these  "interesting''  regions  and  show  how  a 
decomposition  of  CO(M)  can  be  obtained  by  composing  regions  of  local  configuration  spaces, 
CO(Ai,  Aj).  We  begin  by  defining  what  the  region  diagram  of  a  mechanism  is  and  then  show  how 
it  is  possible  to  build  a  region  diagram  by  composing  local  region  diagrams  of  pairs  of  objects.  In 
the  second  section,  we  show  that  for  a  restricted  class  of  mechanisms  (fixed  axis  mechanisms)  the 
configuration  space  can  be  fully  characterized  by  the  cross  product  of  parameterized  set  descrip- 
tors and  a  set  of  equations  relating  their  parameters.  Therefore,  it  is  possible  to  use  a  simplified 
version  of  the  composition  operation  to  put  together  local  regions  into  global  regions.  Finally,  we 
give  an  algorithm  for  the  computation  of  the  regions  of  the  configuration  space  of  this  class  of 
mechanisms. 

Definition  2.0.1:  Two  disjoint  connected  regions  Ki  and  Kj  of  CO(M)  are  contiguous  if  there 
is  a  motion  it  of  M  from  a  point  p0  in  Ki  to  a  point  pf  in  Kj,  that  does  not  pass  through  any  point 
outside  Ki  or  Kj, 

contiguous(Ki,  Kj)  o 

3  p0  €  Ki  and  3  p,  €  K  exists-path(M,  p0,  p.)  and  path(ir(M,  p0,  pf))  C  Ki  U  Kj 

Definition  2.0.2:  Let  CO(M)  be  the  configuration  space  of  M.  Let  KI,  K2,  ...  Kn  be  a 
decomposition  of  CO(M)  into  n  disjoint,  connected  subsets.  A  region  diagram  D(M)  of  CO(M)  is 
an  undirected  graph  such  that  the  nodes  are  the  regions  and  there  is  an  edge  between  two  nodes 
iff  the  two  regions  are  contiguous: 

1.  V  =  {Kir 

2.  E  -  {<Ki,  Kj>  |  Ki  and  Kj  are  contiguous  in  CO(M)} 

The  transition  between  one  res.on  and  the  other  can  be  expressed  as  a  set  of  algebraic  boun- 
dary conditions  on  the  parameters  of  Cspace(M).    These  boundary  conditions  are  associated  with 


2.  We  only  represent  here  the  connected  region  of  the  configuration  space  in  which  the  initial  placement  of  the 
mechanism  is  included. 
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the  edge  in  D(.M)  that  indicates  a  transition  between  the  two  regions.    The  initial  region  is  the  one 
that  contains  the  point  p  .    For  example,  the  region  diagram  of  the  pair  in  Figure  2.1  is: 


Rl     ) 1    R2 

Tl 

where  the  condition  for  a  transition  Tl  from  Rl  to  R2  is: 

I  X,    "    XB  I   =       1, 
I  YA    -    YB  I   =       ° 

2.1.  Local  Region  Diagram 

A  local  region  diagram  of  a  subset  M'  of  M  is  a  region  diagram  for  M'  of  CO(M')  indepen- 
dent of  the  rest  of  the  parts.  In  particular,  we  are  interested  in  local  region  diagrams  of  pairs  of 
objects.  We  consider  first  the  mechanism  formed  only  by  A  and  B,  VT  =  {A,  B}  as  an  indepen- 
dent mechanism  whose  configuration  space  CO(A,  B)  (and  not  COM(A,  B))  is  12  dimensional 
with  respect  to  a  fixed  reference  frame  F. 

2.2.  Region  Classes 

There  are  certain  pairwise  interactions  that  occur  frequently  in  mechanisms  such  as  prismatic 
pairs,  revolute  pairs,  parallel  gears,  etc.  As  we  will  see  below,  each  such  interaction  corresponds 
to  a  region  of  the  pairwise  configuration  space  that  can  be  characterized  by  a  cross  product  of  the 
parameterized  set  descriptors  introduced  in  the  previous  chapter.  Each  pair  defines  a  region  class. 
Following  is  the  description  of  these  classes. 

Let  p  =  <xi,  yi>  ,  1  s  i  <  6  be  a  point  in  CO(A,  B)  and  let  pQ  be  the  initial  placement.  We 
define  the  following  interesting  region  classes  for  pairs  of  objects  that  describe  the  6  lower  pairs 
introduced  by  Reuleaux  [Reuleaux  1876]  (see  Figure  2.2)  and  other  special  higher  pairs  such  as 
gear  pairs. 

Definition  2.2.1:  The  Prism  class 

The  region  R  of  CO(A,  B)  belongs  to  the  prism  class  iff  R  can  be  described  by  the  following 
possible  paths  of  the  reference  points  Oa  and  Ob  with  respect  to  a  fixed  reference  frame. 

Prism 

p-translation(A,  O,  XA)  x  p-translation(B,  O,  XB), 
Lmin  s  |  XA  -    XB  |  ^  Lmax 

where  X  and  X  are  measured  from  the  same  origin  in  O  and  O  is  the  axis  of  translation.  Note 
that  since  p-translation(A,  O,  XA)  is  a  subset  of  COM(A)  and  p-translation(B,  O,  XB)  is  a  subset 
of  CO  (B),  their  cross  product,  plus  the  condition  Lmin  s  |  X^  -  Xg  |  <  Lmax  is  a  subset  of 
CO(A,  B)  (here  M  =  {A,  B}).  In  this  case  R  €  Prism  and  we  write  R  =  prism(A,  B,  O).  The 
interpretation  is  that  A  and  B  can  have  possible  independent  translations. 

We  can  define  in  a  similar  fashion  the  other  lower  pairs  classes:  Revolute,  Cylinder,  Sphere, 
Screw  and  Planar.  Special  pair  classes,  such  as  meshed  parallel  gears  can  also  be  defined: 

Definition  2.2.2:  The  Parallel-Gears  class 

Two  objects,  A  and  B  form  a  parallel  gear  pair  if  their  only  possible  motion  are  two  coordinated 
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rotations: 

Parallel  Gears 

p-rotation(A,  O,  9  J  x   p-rotation(B,  O',  6B), 

9     =   -  r    /  r    x  9  where  r     and  r    are  the  radii  of  A  and  B,  respectively. 

In  this  case,  R  =  parallel-gears(A,  B). 

2.3.    Fixed  Axis  Mechanisms 

A  mechanism  is  a  fixed  axis  mechanism  iff  all  the  possible  motions  of  the  parts  can  be 
described  uj  motions  along  an  axis  that  is  fixed  in  space.  Thus  in  a  fixed  axis  mechanism,  the 
rotary  axis  do  not  move  in  space.  The  piston  assembly  in  Figure  2.4  is  not  a  fixed  axis  mechanism 
since  the  motion  of  the  rod  R  cannot  be  described  as  a  motion  around  or  along  an  axis  that  is 
fixed  in  space.  More  formally,  M  is  a  fixed  axis  mechanism  iff  for  all  its  parts  Ai, 

n 

COM(Ai)  =    (J  simple-region(Ai,  Ok,  XJ 
t-i 

where  a  simple  region  is  a  connected  set  that  can  be  represented  by  one  of  the  following 
parametrized  set  descriptors  corresponding  to  possible  simple  motions: 

1.  c  =  fixed(Ai,  p) 

2.  c  =  p-rotationM(Ai,  0|k,  9|k),  9|k  £  [9min,  9max] 

3.  c  =  p-translationM(Ai,  0|k,  Xfc),  Xfc  €  [Xmin,  Xmax] 

4.  c  =  p-cylinderM(Ai,  Ok,  9|k,  X;k),  X|k  €  [Xmin,  Xmax],  6|k)  €  [6min,  9max] 

5.  c  =  p-helical^Ai,  0(k,  9ik,  X|k),  6ik  =  f(X(k ),  9jk  €  [6min,  9max] 

that  is,  a  region  is  simple  iff  the  only  possible  motions  are  simple  motions  along  the  specified 
axis,  or  no  motion  at  all  (fixed(Ai,  p)).  To  every  parameter  there  is  an  interval  associated  with  it, 
interval(X   ),  possibly  a  function  of  the  other  parameter  as  in  cases  4  and  5. 

Definition  2.3.1:  A  configuration  space  of  a  single  part  Ai,  COM(Ai)  is  a  simple  configura- 
tion space  iff  it  is  the  union  of  n  simple  regions, 


CO    (Ai)  [J     c  where  c    is  a  simple  region 


'  =  i 


Note  that  not  all  unions  of  simple  configuration  spaces  are  possible  (physically  realizable)  confi- 
guration spaces.  Nevertheless,  since  we  start  with  actual  physical  objects,  we  will  end  up  only  with 
feasible  configuration  spaces. 

A  fixed  axis  mechanism  has  a  simple  configuration  space  that  can  be  described  as  a  union  of 
simple  regions  of  the  form: 

•  simple-region(Al,  0lkl,  X|kl)  |x  |  .  .  .  .  |x  |  simple-region(An,  Onkn,  XnkJ 

•  intervals  for  each  parameter,  X    €  [Xmin  ,  [X™  1 

•  A  function  (possibly  discontinuous)  F(XU,  ...  XJnl,  Xnkn)  ^  0 

Definition  2.3.2:  The  decomposition  of  CO(M)  into  a  set  of  simple  regions  is  maximal  iff  no 
two  regions  can  be  merged.  Two  regions  Rl  and  R2  can  be  merged  iff: 
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1.  All  the  parametrized  set  descriptors  specify  the  same  possible  motion,  along  the  same  axes 
and  with  the  same  parameters. 

2.  The  F  function  for  both  regions  is  identical 

3.  The  intersections  between  intervals  are  all  non-empty 

We  want  a  maximal  decomposition  of  CO(M)  since  this  guarantees  that  the  poss.ble  motions  in 
each  region  are  all  qualitatively  different. 

The  advantage  of  simple  regions  is  that  the  composition  of  rwo  simple  regions  yields  always 
a  new  set  of  simple  regions.  Since  the  axis  are  fixed  in  space,  it  is  possible  to  take  the  same  fixed 
reference  frame  to  describe  local  configuration  spaces  and  compose  them.  This  is  not  the  case  in 
movable  axis  mechanisms,  since  we  cannot  use  a  fixed  reference  frame  and  still  be  able  to 
describe  the  local  configuration  spaces  as  simple  regions. 

Note  that  by  the  definition  of  the  fixed  axis  mechanism,  there  is  associated  to  every  object  a 
set  of  axis  along  which  its  configuration  space  can  be  described.  If  we  were  given  in  advance  a  set 
of  fixed  axis  for  every  pair  of  objects,  we  could  compute  four  dimensional  configuration'  spaces 
(one  for  each  axis,  with  rwo  parameters  for  each  object)  instead  of  dealing  with  a  12  dimensional 
configuration  space.  As  the  following  theorem  shows,  given  a  complete  set  of  axis  for  every  pair, 
it  is  possible  to  compute  exactly  the  configuration  space  of  the  mechanism  by  using  a  simplified 
version  of  the  composition  operator. 

Definition  2.3.3:  The  set  of  axis  {O  }  is  a  complete  set  of  axis  with  respect  to  Ai  iff  COM(Ai) 
can  be  completely  characterized  by  a  set  of  simple  regions  along  O  .  Similarly,  the  set  {OJ  is  a 
complete  set  of  axis  with  respect  to  a  subassembly  M'  iff  its  configuration  space  can  be  completely 
characterized  by  a  set  of  simple  regions  along  these  axis. 

Theorem  5:  Let  M  be  a  fixed  axis  mechanism  and  CO(M)  its  corresponding  simple  confi- 
guration space.  Let  {O  }  be  a  given  complete  set  of  fixed  axis  for  every  Ai.  Then  the  computa- 
tion of  CO(M), 

CO(M)  =  CO(Al,  A2)  |x  |  CO(Al,  A3)  |x  | .  .  .  .  |x  |  CO(An-l,  An) 

COM(Al,A2)|x|COM(Al,A3)|x| |x  |  COM(An-l,  An)      (3) 

can  be  carried  using  a  simplified  version  of  the  composition  operation  |x  |  along  the  sets  of  axis. 

Proof:  Let  CO'(Ai,  Aj)  be  the  configuration  space  of  Ai  and  Aj  with  respect  to  a  complete 
set  of  fixed  axis  {OJ.  Since 

COM(Ai,  Aj)  C  CO'(Ai,  Aj)  C  CO(Ai,  Aj) 
and  (3)  is  true  by  Theorem  2  and  the  monotonicity  of  proprerty  of  |x  |  (property  8), 
CO(M)  =  CO'(Al,  A2)    |x|   CO'(Al,  A3)  |x  |   CO'(An-l,An) 

Now  we  show  how  to  compute  the  composition  of  simple  pairwise  configuration  spaces. 

Let  CO'(Ai,  Aj)  =  {c.}  and  CO'(Ai,  Ak)  =  {d^}  be  two  decompositions  of  pairwise  configuration 
spaces  into  simple  regions.    Every  simple  region  is  of  the  form: 

simple-motion(Ai,  O ,,  X^)  X  simple-motion(Aj,  Oy,  X    ). 

XAj  €  [Xmin,  Xmax] 

x'  =  f(X.  ) 

Ai         v    Ar 

where  X     is  the  appropriate  parameter  vector. 

The  composition  of  the  rwo  configuration  spaces,   CO'(Ai,   Aj)   |x  |  CO'(Ai,   Ak)   leads  to  new 
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simple  regions  of  the  form: 

simple-region(Ak,  O,,  X^k)  x  simple-region(Aj,  O,,  Xv) 

x 

(simple-region(Ai,  Or  X,  )  PI  simple-region'(Ai,  O  ,  X^)) 

XAi,  XAj,  XAk  €  Intervals 

*A!  =  f(x,;)'  x„  =  g(x^) 

and  appropriate  boundary  conditions  on  the  parameters  of  possible  motions  derived  from  the  ori- 
ginal intervals  and  the  relations.  The  intersection  rules  for  possible  simple  motions  are  straight- 
forward and  are  listed  in  the  next  section:  the  intersection  of  two  rotations  is  a  rotation  only  if  the 
rotations  is  along  the  same  axis;  the  intersection  between  two  translations  is  not  empty  only  when 
the  two  axis  of  translation  are  parallel,  etc. 

Note  that  the  composition  of  pairwise  simple  configuration  spaces  of  one  parameter  yields  a 
simple  configuration  space  since  the  intersection  of  two  simple  regions  is  a  simple  region,  and  the 
cross  product  of  simple  regions  is  a  simple  region.  Therefore,  to  compute  the  composition  of  sim- 
ple configuration  spaces,  it  is  sufficient  to  be  able  to  compute  the  intersection  of  simple  regions 
and  the  new  intervals  for  each  parameter.  The  intersection  of  two  two-parameter-configuration 
spaces  or  a  two-parameter  and  a  one-parameter  configuration  space  can  yield  a  new  set  of 
regions  which  is  not  connected.  Nevertheless,  the  connected  components  are  simple  regions  (see 
Figure  2.5). 


2.3.1.    Intersections  of  simple  regions 

Let  A  be  a  part  of  M,  O  and  O'  two  axis.  Let  X,  X',  6  and  6'  parameters  and  let  interval(X) 
be  the  interval  associated  with  parameter  X.  Let  relations(X)  ::  X  =  F(Z)  be  the  function  that 
relates  the  parameter  X  of  A  with  the  parameters  of  the  other  parts  (denoted  by  the  parameter 
vector  Z).    The  intersection  of  two  simple  regions  c  and  d  of  COM(A)  is: 

1.  p-rotationM(A,  O,  9)  n  p-rotationM(A,  O,  6')  =  p-rotationM(A,  O,  0") 
interval(9")  =  interval(e)  n  interval(d') 

9"  =  relations(9)  and  9"  =  relations(9') 

2.  p-translationM(A,  O,  X)  n  p-translationM(A,  O',  X')  = 

if  parallel(0,  O')  then  p-translationM(A,  O,  X") 
interval(X")  =  interval(X)  n  interval(X') 
X"  =  relations(X)  and  X"  =  relations(X') 

3.  p-cylinderM(A,  O,  9,  X)    f~l    p-rotationM(A,  O,  9')  = 
p-rotationM(A,  O,  9") 

interval(9")  =  interval(9)  n  interval(9') 
9"  =  relations(9)  and  9"  =  relations(9') 

4.  p-cylinderM(A,  O,  9,  X)    n    p-translationM(A,  O',  X')  = 

if  parallel(0,  O')  then 
p-translationM(A,  O,  X") 
interval(X")  =  interval(X)  D  interval(X') 
X"  =  relations(X)  and  X"  =  relations(X') 

5.  p-cylinderM(A,  O,  9,  X)    f~l    p-cylinderM(A,  O',  9',  X')  = 

if  O  =  O'  then 
p-cylinderM(A,  O,  6",  X") 
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interval(X")  =  interval(X)  n  intervaKX") 
interval(e")  =  interval(9)  fl  interval(8') 
8"  =   relations(8)  and  8"  =   relations(9') 
X"  =  relations(X)  and  X"  =  relations(X') 

if  paralleKO,  O')  then 
p-translationM(A,  O,  X") 
interval(X')  =  interval(X)  f~l  interval(X') 
X"  =  relations(X)  and  X"  =  relations(X') 

6.  p-helicalM(A,  O,  8,  X)  O  p-hehcalM(A,  O,  8',  X')    = 

if  8  =  8'  then 

p-helicaL^A,  O,  9",  X") 

interval(8")  =  interval(8)  D  interval(8') 

7.  In  all  other  cases,  or  when  the  intersection  of  parameter  intervals  is  empty  or  the  relations 
are  incompatible,  the  intersection  is  one  single  point,  corresponding  to  the  object  A  being 
fixed  in  space. 

The  intersection  of  intervals  for  the  simple  motions  p-rotation  and  p-translation  is  straight- 
forward since  it  consists  of  intersecting  two  linear  intervals  and  checking  the  consistency  of  the 
corresponding  equations.  The  intersection  of  two  intervals  of  angles  can  yield  at  most  two  regions 
and  can  also  be  computed  easily.  This  is  not  the  case  however  for  p-cylindrical  since  the  boun- 
daries can  be  expressed  as  a  function  X  =  f(8)  instead  of  as  a  simple  interval.  In  this  case  we 
need  to  intersect  functions  to  derive  a  new  function  that  describes  the  bounds.  Section  2.6  will  dis- 
cuss how  to  carry  this  operation  and  how  to  characterize  the  disconnected  pieces. 

2.4.    Computing  the  region  diagram  of  fixed  axis  mechanisms 

Theorem  5  suggests  a  direct  method  of  computing  CO(M)  for  fixed  axis  mechanisms:  for 
each  pair  of  objects,  select  a  set  of  potential  axis,  compute  their  configuration  space  with  respect 
to  these  axes  and  then  compose  the  pairwise  configuration  spaces  to  obtain  CO(M). 

It  is,  however,  not  necessary  to  first  compute  all  pairwise  configuration  spaces  and  then  com- 
pose them.  The  computation  of  certain  pairwise  configuration  spaces  is  simplified  if  we  know 
what  are  the  possible  motions  of  the  objects.  For  example,  if  it  is  known  that  A  is  fixed  and  that 
B  can  only  translate  along  an  axis  O,  then  if  the  motion  envelope  of  B  does  not  intersect  A, 
CO'(A,  B)  can  be  simply  ignored  in  the  composition.  It  therefore  makes  sense  to  attempt  to  com- 
pute some  pairwise  configuration  spaces  first,  and  then  use  the  information  derived  from  them  to 
simplify  (or  even  ignore)  the  computation  of  other  pairwise  configuration  spaces. 

The  above  remarks  suggest  the  following  algorithm  to  compute  CO(M):  First,  compute  the 
contact  graph.  In  this  graph,  every  object  is  represented  by  a  node  and  there  is  an  edge  between 
two  nodes  iff  the  two  corresponding  objects  are  in  contact.  Then,  find  the  set  of  objects  in  contact 
with  the  fixed  supporting  frame  and  find  the  possible  motions  of  these  objects  by  analyzing  the 
pair  (frame,  object).  Label  every  node  of  the  graph  with  the  possible  motion  found  for  the 
corresponding  object.  Next,  find  the  set  of  objects  in  contact  with  the  ones  already  labeled  and  do 
a  pairwise  analysis,  adjusting  and  adding  the  possible  motions  accordingly.  Continue  doing  the 
pairwise  analysis  of  objects  in  contact  until  no  more  changes  are  made  to  the  labels.  Using  the 
labels  in  each  node,  build  a  portion  of  the  region  diagram.  Finally,  determine,  using  the  labels  of 
each  node,  if  new  contacts  can  be  created  or  old  contacts  can  be  interrupted.  If  so,  update  the 
contact  graph,  analyze  the  new  pairs  and  update  the  labels.  Repeat  this  process  until  all  possible 
topologies  are  explored.  At  the  end,  the  region  diagram  will  be  constructed.  This  method  is  simi- 
lar to  doing  constraint  propagation,  where  the  possible  labels  are  the  set  descriptors  p-rotation,  p- 
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translation,  etc.  and  the  intervals  of  the  parameters,  and  the  constraints  are  the  equations  relating 
the  parameters.  The  constraint  propagation  is  done  for  each  contact  topology,  and  the  region 
diagram  is  then  extracted  fom  it. 

To  describe  the  algorithm  in  detail,  we  introduce  the  following  notation: 

•  C(M)  =   G(V,  E)  is  the  contact  graph.  The  function  neighbors(V)  returns  the  set  of  nodes 
connected  to  node  V. 

•  F  is  the  fixed  supporting  frame  of  the  mechanism  M. 

2.4.1.  Algorithm  1:  Sketch 

1.  Build  a  contact  graph  C(M)  for  the  mechanism  M  in  the  initial  placement. 

2.  A  :=  F  /*  A  is  the  first  object  in  the  pair)  */ 

S  :=  neighbors(F)     /*  S  is  the  set  of  objects  in  contact  with  F  */ 
S'  :=  0  /'  S'  is  the  set  of  nodes  left  to  be  updated  V 

3.  DO 

3.1  For  all  B  in  S  do 

3.1.1  Find  a  set  of  possible  fixed  axis  for  the  pair  (A,  B) 

3.1.2  Compute  CO'(A,  B),  taking  into  account  possible  fixed  axis 
and  the  labels  of  A  and  B. 

3.1.3  Update  the  labels  of  A  and  B  according  to  CO'(A,  B). 
If  either  changed,  add  it,  and  its  neighbors  to  the  set  of 
nodes  to  be  updated. 

S'  :=  S'  U  neighbors(B)  U  neighbors(A) 

3.2  End  do 

3.3  A  :=  first(S')  /*  Choose  a  new  first  object  for  the  pair  */ 

S  :=  S'  -  A  /•  Update  S'  */ 

Until  S'  is  empty 

4.  Construct  the  part  of  the  region  diagram  that  corresponds  to  the  labeling 
in  the  graph. 

5.  For  every  pair  (Ai,  Aj)  not  in  contact  in  C(M),  check  whether  they  can 
come  in  contact  given  the  labels  on  Ai  and  Aj.  If  they  can,  modify  the 
contact  graph  appropriately. 

6.  For  every  pair  in  contact,  check  whether  they  can  stop  being  in  contact 
and  modify  the  contact  graph  appropriately. 

7.  Repeat  steps  2  to  6  until  all  possible  topological  changes  have  been  explored. 

Claim:  The  above  algorithm  is  complete,  i.e  it  will  compute  CO(M)  correctly  iff  the  set  of 
potential  axis  for  each  pair  includes  a  complete  set  of  axis  for  the  pair. 

2.4.2.  Finding  Possible  Fixed  Axis 

To  find  possible  candidates  for  fixed  axis,  we  note  that  by  the  definition  of  lower  pairs,  two 
objects  form  a  lower  pair  iff  they  have  a  surface  contact.  Therefore  natural  candidates  for  axes 
are: 


1. 
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For  planar  surfaces,  axes  parallel  to  the  surface  of  contact  and  to  the  edges  in  contact  with 
surfaces. 


2.       For  cylindrical  surfaces,  axes  that  coincide  with  the  axis  of  the  cylinder 

For  each  pair  of  objects,  the  set  of  possible  axes  can  easily  be  derived  from  an  initial  placement  of 
objects  in  contact  and  the  geometrical  description  of  each  object.  Therefore,  we  can  guarantee  that 
COM(A,  B)  C  CO'(A,  B). 

2.4.3.    Example:  Region  Diagram  for  the  Clutch  System 

In  this  section  we  will  illustrate  how  to  build  the  region  diagram  of  the  clutch  system  in  Fig- 
ure 2.5.  This  system  works  as  follows:  F  is  the  fixed  supporting  frame,  A  is  the  input  shaft  and  C 
is  the  output  shaft.  Initially  A  and  C  are  not  in  contact,  so  that  the  rotation  of  A  is  not  transmitted 
to  any  part.  Mounted  on  the  output  shaft  C  is  a  disk  B  which  is  mounted  on  C  on  a  square  cross 
section  (see  the  projections  in  the  bottom  of  Figure  2.5).  This  means  that  B  can  slide  back  and 
forth,  and  that  when  B  turns,  C  also  turns.  Mounted  also  on  C  is  part  D,  which  slides  together 
with  B  but  does  m  rotate.  The  clutch  is  engaged  by  screwing  E,  which  causes  D  to  slide,  pushing 
in  turn  B.  If  the  4  cylinders  mounted  on  B  are  positioned  to  match  the  holes  in  the  disk  of  A,  then 
A  and  B  fit  together.  Thus,  a  rotation  of  A  will  cause  a  rotation  of  B,  which  in  turn  will  make  the 
output  shaft  rotate  by  the  same  speed  and  in  the  same  direction.  The  clutch  is  disengaged  by 
unscrewing  E. 

Following  Algorithm  1,  the  analysis  of  the  clutch  goes  as  follows:  we  first  construct  the  con- 
tact graph: 


Let  F  be  the  fixed  reference  frame  of  the  mechanism.  We  start  by  examining  all  the  neighbors  of 
F:  A,  C,  D,  E.  For  each  of  the  pairs  <F,  A>,  <F,  C>,  <F,  D> ,  and  <F,  E>  we  analyze  their 
configuration  space.  Since  all  are  cylinders,  their  potential  natural  axes  are  their  central  axes. 
Thus,  the  analysis  of  their  configuration  spaces  along  the  fixed  axis  O  (for  A  and  C)  and  of  O' 
(for  E)  will  yield  that  both  <F,  A>  and  <F,  C>  form  a  revolute  pair  and  that  <F,  E>  forms  a 
helical  pair.  Since  the  frame  F  is  fixed,  the  possible  motions  are: 

fixed(F,  p) 

p-rotation(A,  O,  9A) 

p-rotation(C,  O,  8C) 

p-helical(E,  O',  X£,  6E),    X£  <E  [0,  Xmax],  X£  =  pitch  *  9£ 

As  for  the  analysis  of  the  pair  <F,  D> ,  there  are  several  potential  axes:  the  axis  of  the  cylinder, 
O,  and  an  axis  perpendicular  to  O,  since  there  is  a  surface  contact  between  F  and  D  in  that  direc- 
tion. Although  F  and  D  form  a  plannar  pair,  the  algorithm  can  only  analyze  pairs  along  specified 


0. 
1. 
2. 
3. 


6,  i  [0,2^] 
6C  €  [0,  2ir] 
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axes.  Thus  it  will  analyze  F  and  D  along  two  candidate  axis,  O  and  O".  The  analysis  along  these 
two  axes  yields  that  the  pair  <F,  D>  is  a  prismatic  pair  along  both  axes. 

4.  p-translation(D,  O,  Xc),     X£  €  [0,  Xmax] 

5.  p-translation(D,  O",  X'E),  X'E  t  [0,  X'max] 

We  now  look  at  objects  in  contact  with  A,  C,  D,  E  other  than  F.  We  find  the  pairs:  <E,  D> , 
<D,  C>,  <D,  B>  and  <B,  C>  The  analysis  of  the  pair  <E,  D>  is  done  along  the  potential 
axes  previously  found:  O,  O'  and  O"  and  taking  into  account  the  existing  labels  for  E  and  D.  The 
intersection  of  3  and  4  and  of  3  and  5  yield  the  following  new  labels: 

3.        p-helical(E,  0\  X£,  8E),    XE  €  [0,  Xmax],  X£  =  pitch  *  % 

6.  p-translation(D.  O,  XD),  XD  €  [0,  Xmax],  XD  =  X£ 

Continuing  this  process  we  obtain  the  following  complete  labeling: 

fixed(F,  p) 

p-rotation(A,  O,  9A) 

p-rotarion(C,  O,  9C) 

p-helical(E,  0\  X..  8_), 

XE   =  p.tch(E)  -9CE,X^  <  XEN0AGE 

p-translation(D,  O,  XD), 

xD  =  xE 


p-cylinder(B,  O,  Xg,  6B) 

Xg  =  xD 

6r 


B  C 

This  corresponds  to  a  single  region,  DISENGAGED. 

The  next  step  is  to  check  if  there  can  be  new  contacts  created.  By  examining  the  pairwise 
possible  motions  and  their  geometry,  we  determine  that  only  A  and  B  can  come  into  contact.  See 
Figure  2.6  for  the  configuration  space  of  A  and  B  along  axis  O.  There  are  4  new  regions 
described  as: 


ions  2  to  5:  ENGAGED  i 

fixed(F,  p) 

p-rotation(A,  O 

•9a) 

p-rotation(C,  O, 

ec) 

eA  =  ec  *  i*iT/: 

2,            1  £ 

i  s  4 

p-helical(E,  O', 

XE,8E), 

X^  =  pitch(E)  ' 

9E,  X£  > 

Y 

ENGAGE 

p-translation(D, 

o,xD), 

XD   =    X  E 

p-cylinder(B,  O. 

XB,  9B), 

XB=    XD 

6B=    6C 

The  transitions  between  the  regions  are: 


Transitions 


eA  =  9C  -  i>/2 

XB   >    XENGAGE 
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The  complete  region  diagram  is  thus: 


^DISENGAGE 


This  corresponds  to  the  qualitative  description  of  the  possible  behaviors  of  the  clutch. 

2.5.    Regions  within  Regions 

Let  A,  B  be  a  pair  of  objects  and  let  R  be  a  simple  connected  region  of  their  configuration 
space  CO(A,  B).  Since  R  is  a  simple  region,  its  most  general  form  is: 

simple-motion(A,  Ol,  Xv  9^)  x  simple-motion(B,  02,  Xg,  9B) 

f(xv9A,xB,eB)^o 


XA  €  [XminA,  Xm"A] 
XB  €  [X«"*"B,  Xm"B] 

eA  €  [9T,nA>  emaxA] 
eR  €  (eminR,  emax„] 


'B         L  B' 

where  f  is  a  function,  possibly  discontinuous  on  the  intervals  of  the  parameters.  Although  R  is  a 
simple  connected  region,  it  can  in  fact  reflect  two  qualitatively  different  behaviors.  This  will  be 
illustrated  in  the  following  example. 

Let  A  be  a  rod  which  is  free  to  rotate  around  Ol,  and  let  B  be  a  square  that  is  free  to 
translate  along  axis  02,  as  shown  in  Figure  2.7.  The  region  diagram  consists  of  two  disconnected 
simple  regions.  The  upper  one  is  described  by: 

p-rotation(A,  Ol,  6A)    |x  j  p-translarion(B,  02,  Xn,  Xn) 

8A  €  [0,  2tt],  Xb  €  [0,  lmax] 

f(8    ,  X  )  s  0,  where  f  is  as  shown  in  the  diagram 

Now  suppose  that  the  block  B  is  always  kept  in  contact  with  A.  A  rotation  of  A  from  0  to  it 
implies  a  translation  of  B  upwards  and  then  a  translation  of  B  downwards.  Thus,  there  are  two 
qualitatively  different  behaviors  that  can  be  achieved  in  this  region.  Since  our  goal  is  to  have 
every  region  in  the  region  diagram  reflect  one  simple  type  of  qualitative  behavior,  we  need  to 
further  partition  the  region  into  subregions  that  will  reflect  qualitative  changes  of  possible 
behavior.  In  the  previous  case  we  can  partition  the  region  into  the  following  three  regions: 

Rl:  8A  €  [0,  tt/2],  Xr  >  fl(6A),  XB  s  lmax 
R2:  8^  €  [tt/2,  tt],  Xs  ^  f2(9J,  Xn  <  lmax 
R3:    9^  €  [tt,  2tt],  Xb  >  0,  X3  <  lmax 
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where  fl  is  monotonically  decreasing  and  f2  is  monotonicaliy  decreasing.  Note  that  this  is  one 
possible  decomposition  of  the  region  along  the  9A  parameter.  A  different  decomposition  can 
result  if  we  decompose  along  the  Xg  parameter  (the  6A  decomposition  is  more  useful  if  the  rod  is 
driving  the  block  and  the  X3  decomposition  is  more  useful  if  the  opposite  is  true). 

Definition  2.5.1:  Decomposition  of  Simple  Regions 

A  decomposition  of  R  into  subregions  is  a  set  C  =  {C{,  ...  C^}  such  that: 

1.  All  the  subregions  Ct  are  pairwise  disjoint. 

2.  The  union  of  all  subregions  C  's  is  R. 

3.  Let  f(X   ,  6    ,  XD,  9D)  2  0  be  the  function  that  defines  the  region  C  and  let  f  (X    ,  9    ,  X,, 

v  A  A  o  a  I  A  A3 

9  )>  0  be  the  piece  of  f  that  defines  subregion  C  .  Then,  for  every  subregion,  the  function 
f  (X  ,9  ,  X„,  93)  a=  0  can  be  decomposed  into  a  set  of  functions  f_  such  that  every  such 
function  is  both  continuous  and  monotonic  over  the  intervals  of  the  subregion. 

4.  The  decomposition  is  maximal,  i.e.  no  two  subregions  can  be  merged  into  another  subregion 
satisfying  properties  3  and  4. 

The  purpose  of  condition  3  is  to  guarantee  that  the  boundaries  of  a  subregion  vary  only  in 
one  direction  (either  monotonically  decreasing,  increasing  or  constant")  so  that  this  will  correspond 
to  a  monotonic  movement  if  the  movement  is  a  ong  the  boundary.  Note  that  the  definition  allows 
several  possible  decompositions  of  a  region.  For  the  previous  example,  Figure  2.8  shows  two  pos- 
sible decompositions  of  R:  one  along  9  consisting  of  the  boundary  portions  fn,  f,,  and  f31  and  one 
along  X  whose  boundary  portions  are  fM,  fp,  f:],  f„,  and  fjr  In  general,  for  every  region  R  there 
are  four  possible  (not  necessarily  distinct)  decompositions  of  R. 

Definition  2.5.2:  Let  interval(X,  C  )  be  the  interval  of  parameter  X  in  subregion  C  .    We  say 
that  a  decomposition  C  of  R  is  along  parameter  x  iff  all  the  intervals  are  pairwise  disjoint,  i.e 
interval(X,  C)  n    interval(X,  C)  =  0,     i  *  j 
and 

n 

1J     interval(X,  C  )  =  interval(X,  R)      for  all  i 

i=  i 

We  extend  the  concept  of  subregions  to  regions  of  the  full  configuration  space  CO(M)  by 
extending  the  requirements  over  the  four  parameter  function  f  to  the  the  n  parameter  function 
describing  the  region. 
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3.    State  Diagram 

Given  a  region  diagram  for  a  mechanism,  its  initial  placement  and  a  sequence  of  input 
motions  applied  to  its  various  parts,  we  want  to  find  what  are  the  actual  motions  of  all  the  parts  in 
the  mechanism  and  their  relationships.  Furthermore,  we  want  a  qualitative  distinction  between  the 
different  motions  that  each  part  has  and  the  different  states  through  which  the  mechanism  goes 
through.  Our  notion  of  object  states  is  similar  to  the  one  presented  in  [Forbus  85]  and  [Hopcroft 
and  Krafft  87]:  the  states  of  an  object  correspond  to  qualitatively  different  behaviors.  The  connec- 
tions and  transitions  between  the  states  of  each  object  can  be  described  in  a  state  diagram.  The 
basic  distinction  between  the  state  diagram  and  the  region  diagram  is  that  the  region  diagram 
deals  with  possible  motions  of  parts,  while  the  state  diagram  deals  with  actual  motions  due  to 
some  input  motion  applied  to  parts  of  the  mechanism.  The  state  diagram  reflects  the  qualitative 
behavior  of  the  mechanism  for  a  given  input  sequence. 

3.1.    Input  Motions 

We  restrict  the  possible  input  motions  to  be  a  sequence  of  simple  motions:  rotation,  transla- 
tion and  helical  motion  around  a  fixed  axis.  Formally,  let  T  be  the  input  motion  alphabet  such  that 
r  =  {rotation(A,  O,  9),  translation(A,  O,  X),  helical(A,  O,  X,  9)}  where  A  is  a  part  of  M,  O  is  a 
fixed  axis  and  each  parameter  X  has  specified  with  it  an  interval  X  €  [Xmin,  Xmax],  and  a  direc- 
tion of  change,  direction(X)  =  -  or  -  ,  depending  on  whether  X  varies  from  Xmin  to  Xmax  or 
the  opposite.  In  the  case  of  helical  motion,  we  also  specify  the  relationship  between  9  and  X  by  a 
monotonic  function  f  over  the  given  interval  for  X,  X  =  f(6). 

Definition  3.1.1:  A    simple  input  motion  sequence  is  a  sequence  of  simple  motions  a. 

For  example,  an  alternate  translation  of  A  is  represented  by: 

cr  =  [translation(A,  O,  xl),  translation(A,  O,  x2)] 

where 

xl,  x2  €  [Lmin,  Lmax]  and  direction(xl)  =   +,  direction(x2)  =   - 

Note  that  this  definition  implies  that  the  motions  are  applied  one  after  the  other  so  that  no  two 
motions  can  be  applied  in  parallel  to  two  different  objects. 

Theorem  5:  Let  M  be  a  fixed  axis  mechanism.  Then  if  the  input  motions  form  a  simple 
sequence,  the  actual  motions  of  the  parts  of  M  can  also  be  described  by  sequences  of  simple 
motions. 

Proof:  The  argument  is  based  on  the  fact  that  in  a  fixed  axis  mechanism  only  simple  motions 
can  occur,  and  therefore  the  parts  will  either  have  a  simple  motion,  or  no  motion  at  all.  The  fol- 
lowing section  provides  a  set  of  rules  for  the  propagation  of  input  motions  that  generate  only  sim- 
ple motions. 

Note:  We  are  assuming  that  the  motion  sequence  is  a  legal  sequence,  i.e.  that  the  stated  motions 
are  indeed  possible.  For  example,  if  a  part  is  fixed,  then  applying  to  it  a  translation  is  an  illegal 
motion.  This  restriction  is  made  only  to  facilitate  the  presentation  and  simplify  the  arguments.  We 
will  show  that  we  can  do  without  this  restriction  by  saying  that  the  effect  of  applying  an  illegal 
motion  to  an  object  results  in  no  motion  of  the  object.  Similarly,  if  the  interval  of  a  motion  is  not 
within  the  limits  of  a  legal  interval,  we  turn  this  illegal  motion  into  a  legal  one  by  intersecting  the 
two  intervals.  For  example,  if  the  input  motion  is  rotation(A,  O,  x)  and  x  €  [0,  5],  direction(x)  = 
-  and  the  possible  motion  of  A  is  p-rotation(A,  O,  x),  x  €  [2,  10]  then  the  corresponding  legal 
input  motion  is  rotation(A,  O,  x),  x  €  [2,  5],  direction(x)  =  + .  So  for  every  sequence  a  there 
exist  a  sequence  a'  (possibly  empty)  that  is  legal. 


3.2.  State  Diagram 

Definition  3.2.1:  Let  M  be  a  mechanism,  and  a  a  legal  simple  input  motion  sequence.  A 
state  diagram  S(M,  cr)  is  a  directed  graph  G(V,  E)  such  that: 

1.  The  nodes  of  the  graph  correspond  to  a  description  of  the  actual  motion  of  each  part  and 
their  relationships.  States  are  qualitatively  different  in  that  no  two  states  can  be  merged 
(using  a  criteria  similar  to  the  merging  of  regions)  and  the  states  are  maximal. 

2.  The  directed  edges  between  states  carry  the  conditions  on  the  different  parameters  for  the 
transition  to  be  true. 

Note  that  this  definition  does  not  imply  that  there  is  one  state  per  input  motion.  There  can 
be  cases  in  which  an  input  motion  causes  two  qualitatively  different  behaviors  and  thus  produces 
two  states.  On  the  other  hand,  two  input  motions  correspond  to  at  least  two  states  since  by  defini- 
tion the  two  input  motions  are  qualitatively  different.  So  there  are  at  least  as  many  states  as  there 
are  input  motions  in  the  given  sequence. 

Conjecture:  Given  fixed  axis  mechanism  M  and  an  input  motion  ct,  the  corresponding  state 
diagram  is  unique,  provided  that  no  two  states  can  be  merged. 

Since  the  input  motions  are  one-dimensional,  the  resulting  motions  of  each  part  are  also  one- 
Jimensional.  Since  the  difference  between  states  is  based  on  the  direction  of  change  of  each 
parameter,  there  is  only  one  way  of  partitioning  each  resulting  motion.  This  means  that  although 
the  region  diagram  is  not  unique,  we  will  obtain  the  same  state  diagram  for  for  each  region 
diagram  on  a  given  input  sequence. 

The  region  diagram  describes  the  set  of  all  possible  motions  and  their  relationships  that  parts 
in  a  mechanism  can  have.  Therefore,  all  the  possible  behaviors  of  the  mechanism  are  represented 
there.  The  state  diagram,  on  the  other  hand,  describes  the  actual  behavior  of  the  mechanism  for  a 
specific  input  motion.  We  can  thus  use  the  region  diagram  to  deduce  the  state  diagram  for  a  given 
input  motion  by  turning  possible  motions  into  actual  motions,  as  we  will  see  in  the  next  section. 

3.3.  Algorithm  to  construct  the  state  diagram 

The  state  diagram  of  a  fixed  axis  mechanism  M  can  be  constructed  from  the  following: 

1.  Region  diagram  of  M 

2.  Initial  placement  of  M,  p0 

3.  Simple  input  motion  sequence 

The  algorithm  works  by  propagating  an  input  motion  throughout  the  mechanism  and  "inter- 
secting" the  actual  motion  with  the  possible  motions  inside  a  region.  It  can  be  that  this  actual 
morion  requires  a  transition  to  a  new  region,  in  which  case  the  propagation  continues  in  the  new 
region.  The  transmission  of  motion  from  one  part  to  another  occurs  iff  there  is  a  parameter  rela- 
tion between  the  two  possible  motions  of  the  objects.  For  example,  assume  that  A  and  B  are  two 
meshed  parallel  gears.  Then  their  configuration  space  is: 

p-rotation(A,  O,  9A)  x   p-rotation(B,  O',  oB), 

eA,eB  €  [0,2*] 

9A  =  -  RA/RBx  eB 

Suppose  that  the  input  motion  is  rotarion(A,  O,  9X),  9A  €  [0,  2ir  ]  and  direction(9 K)  =  +.  Then 
the  intersection  of  this  motion  with  the  possible  motion  of  A  yields  rotarion(A,  O,  9^),  9A  €  [0, 
2tt]  and  direction(9  )  =  -.  Since  there  is  a  relation  between  a  parameter  of  A  and  a  parameter 
of  B,  the  motion  is  propagated  and  p-rotation(B,  O',  9B)  becomes  rotation(B,  O',  9g),  9g  €  [0, 
2tt]  and  direction(9^)  =   -  .    The  relation  between  the  motions  is  9^  =   -  RA  /  RB  x  9B.    If  there 
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is  no  relation,  then  the  parts  which  have  not  been  assigned  a  motion  will  be  assigned  the  label  no- 
rnotion(A).  This  reflects  the  default  assumption  that  parts  Jo  not  move  unless  motion  is  transmit- 
ted to  them. 

The  algorithm  works  as  follows:  from  the  initial  placement  of  the  mechanism,  find  the  initial 
region  in  the  region  diagram.  Take  the  first  motion  in  the  input  sequence  and  propagate  ;t  inside 
the  region  by  intersecting  it  with  the  possible  motions  of  each  object.  If  all  the  parameter  intervals 
are  within  the  bounds  specified  by  the  region,  create  a  state  and  proceed  to  the  next  input  motion. 
Else,  find  the  bounds  for  the  intervals  that  are  contained  in  the  region,  and  create  a  state.  Find 
which  transition  condition  in  the  region  diagram  is  satisfied  and  go  to  that  region.  Propagate  there 
the  "remaining"  of  the  input  motion  and  create  a  new  state.  This  process  is  repeated  until  all  the 
parameter  intervals  are  within  bounds,  and  for  all  the  motions  in  the  input  sequence.  At  the  end, 
merge  as  many  states  as  possible. 

In  order  to  present  the  algorithm  in  more  detail,  we  introduce  some  notation: 

•  o-  =  To- ,  a-  1,  the  parameter  in  a   is  x  and  its  interval  is  interval(x). 

L      1»  »       nJ>  r  ,  , 

•  The  region  diagram  is  a  graph  with  regions  Rc,  ...  R^.  Given  a  condition  C,  the  function 
successor(R  ,  C)  returns  the  next  region  by  matching  C  with  the  labels  of  the  edges  that  con- 
nect the  region  R   with  other  regions  of  the  region  diagram. 

•  p-motion(A.  R  ,  x)  is  the  possible  motion  of  A  in  region  R  and  let  interval(x)  be  the  inter- 
val range  of  x. 

•  The  function  propagate(o  ,  R  )  returns  the  actual  motion  of  each  part  in  the  mechanism,  and 
the  updated  a  .  If  the  intervals  of  the  actual  motion  is  contained  in  the  possible  motion  then 
a  is  updated  to  empty.  Otherwise,  the  a  interval  is  updated  to  be  interval(x)  -  interval(x), 
upper  or  lower  portion  depending  on  the  direction  of  x 

Algorithm  3:  State  Diagram 

1.  Find  the  initial  region  R0  from  the  initial  placement  pQ. 
Let  R  :=  RQ 

2.  While  ct  is  not  empty,  do 

2.1  s  :=   first(cr);  o-  :=  a  -   s 

2.2  While  s  is  not  empty,  do 

2.2.1  propagate(s,  R) 

2.2.2  Create  a  new  state  and  instantiate  it  with  the  results 
of  the  propagation. 

Create  a  directed  edge  between  the  actual  state  and  the 

new  state. 

The  label  of  that  edge,  C,    is  the  limit  constraint  on  s. 

2.2.3  If  s  is  not  empty,  R  :=    successor(R,  C). 

2.2.4  End 

2.3  Merge  as  many  states  as  possible 

2.4  End 

Algorithm  4:  Propagation  of  motion 

Given  a  simple  input  motion  s  and  a  region  R,  find  the  actual  motions  of  every  object. 

Let  A  be  the  object  referred  to  in  s.  Let  p-motion(A,  x)  be  its  possible  motion  in  R.  Inter- 
sect s  and  p-motion(A,  x)  using  intersection  rules  similar  to  the  ones  used  for  intersecting  simple 
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regions.  If  the  intersection  is  empty,  label  all  actual  motions  as  "no-motion"  (the  case  of  an  illegal 
input  motion).  Otherwise,  taKe  the  difference  of  intervals  and  return  s  updated  with  this  new 
inte-val. 

The  propagation  of  motion  takes  place  through  the  relationships  between  parameters.  This 
process  is  continued  until  all  relationships  have  been  taken  into  account.  Then  all  parts  which  still 
have  p-morion  labels  are  labeled  as  no-mo:ion  (no  motion  was  propagated  to  them  and  thus  the 
default  is  that  they  will  not  move).  Since  the  boundaries  of  all  the  parameters  are  already 
adjusted  to  reflect  possible  motions,  the  propagation  of  morion  needs  only  to  be  in  the  forward 
direction,  with  no  backtracking. 

3.4.    Example:  The  Clutch  System 

Consider  the  region  diagram  of  the  Clutch  system  described  in  the  previous  section.  Assume 
we  are  given  the  following  input  sequence  and  initial  position: 

•  Initial  state:  Clutch  Disengaged 

xB  =  xD  =  xE  =  0 

K    =    63   =    9C    =    9E    =    ° 

•  Input  motion  sequence: 

cr  =  [rotation(C,  O,  9C),  helical(E.  O',  X_,  ec),  rotation(A,  O,  9  J] 

6C  €  [-it/2,  0],  direction(9c)  =    - 

XE  €  [0,  XENGAGE],  direction(XE)  =   ^,X£  =  pitch  *  9E 

9     €  [0,  2"ir/2],  direction(6A)  =    -,  (this  corresponds  to  aligning  A  and  B,  then  screwing  E 

to  engage  A  and  B  and  then  rotating  A  to  transmit  its  motion  to  C). 

We  build  the  state  diagram  following  the  algorithm  presented.  First,  we  determine  that  the 
initial  region  is  Disengaged.  We  start  by  intersecting  rotarion(C,  O,  9C)  with  p-rotation(C,  O,  6C) 
to  obtain  rotation(C,  O,  6C).  Since  9B  =  9C,  we  have  rotation(B,  O,  6g),  which  intersected  with 
p-rotarion(B,  O,  9B)  gives  rotation(B,  O,  0B).  Since    9B  is  not  related  to  any  other  parameter,  the 

propagation  stops.  All  other  parts  are  labeled  as  no-motion(P)3. 

Now  given  the  next  input  motion,  we  determine  that  there  is  no  region  change  since  the  helical 

morion  of  E  is  possible  inside  the  region.  The  propagation  is  done  in  a  similar  way  as  before.  The 

last  motion  requires  a  change  of  state  since  XB  >    XENGAGED-   We  obtain   the   following  state 

diagram: 

State  Diagram 

[ALIGNMENT] 

rotation(C,  O,  8C),  8C  <E  [-ir/2,  0],  direction(9c)  =   - 
rotation(B,  O,  9B),  9B  €  [-ir/2,  0],  direction(9B)  =   - 

ec)  =    9B 
T(ALIGNMENT,  ENGAGEMENT) 

9C   =    -it/2 

[ENGAGEMENT] 

helical(E,  O",  X£,  9£),  XE  6  [0,  XENGAGE],  direction(XE)  =   -  ,  XE  =  pitch  *  9£ 
translation^,  O,  XD),  XD  €  [0,  XENGAGE)],  direcrion(XD)  =   - 


3  For  compactness,  whenever  the  motion  ot"  a  part  P  is  not  stated  explicitly,  it  is  no-motion(P) 
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translation^,  O,  XB),  X3  $  [0,  X. 

XB  =   X  D 
T(ENGAGEMENT,  COUPLED) 
x    =  X     =  X 

^B  D  ENGAGE 

[COUPLED] 

rotation(A,  O,  6A),  0^  <E  [0,  2tt] 
rotation(B,  O,  9B),  6B  £  [0,  2tt] 
rotation(C,  O,  8C),  6C  €  [0,  2-] 


N"C  \GE 


J,  directionfX. 


e 


e. 
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4.  Movable  Axis  Mechanisms 

Movable  axis  mechanisms  have  several  properties  that  make  their  analysis  hard.  First,  the 
composition  of  simple  local  regions  is  not  a  simple  region.  In  the  piston  example  of  Figure  2.3, 
both  the  rod/cylinder  pair  and  the  rod/bar  pair  are  simple  revolute  pairs  and  thus  have  simple 
relative  motions.  Nonetheless,  we  cannot  compose  their  pairwise  cont: juration  space  since  they  do 
not  have  a  common  reference  frame.  The  actual  motion  of  the  rod  is  not  a  simple  motion.  This 
means  that  we  can  no  longer  use  the  algorithm  for  computing  global  regions  presented  in  section 
2.  Also,  we  can  no  longer  consider  pairwise  four  dimensional  configuration  spaces,  but  rather 
have  to  consider  the  full  12  dimensional  configuration  space.  Still,  the  basic  approach  for  building 
the  configuration  space  from  pairwise  local  configuration  spaces  is  valid  for  any  mechanism  since 
all  the  results  in  Chapter  1  hold  with  no  assumptions.  The  difficulty  is  in  composing  them. 

Another  complication  that  arises  is  the  nature  of  the  possible  motions  parts  can  have.  For 
fixed  axis  mechanisms,  we  were  able  to  identify  five  types  of  qualitatively  different  motions  that 
covered  the  whole  range  of  possible  motions.  This  is  no  longer  true  in  movable  axis  mechanisms. 
[Hopcroft  et  al.  84]  showed  that  it  is  possible  to  build  a  planar  linkage  mechanism  that  represents 
a  given  polynomial  P.  A  planar  linkage  mechanism  is  a  mechanism  consisting  of  n  rods  that  can 
be  joined  together  at  their  endpoints  by  pins  about  which  they  can  rotate.  The  pins- and  rods  are 
allowed  to  move  freely  on  the  plane.  By  selecting  an  appropriate  length  for  the  rods  and  a  partic- 
ular set  of  interconnections,  it  is  possible  to  have  the  endpoint  of  a  rod  R  in  the  linkage  describe 
the  curve  defined  by  the  multivariate  polynomial  P(xr...,  xm)  where  x  represents  the  angular 
position  of  link  i.  This  means  that  the  path  described  by  the  reference  point  of  R  can  be  an  arbi- 
trary polynomial.  Therefore,  if  we  consider  only  the  nature  of  this  curve  as  we  did  for  fixed  axis 
mechanisms,  we  end  up  with  infinitely  many  qualitatively  different  paths  and  thus  infinitely  many 
qualitatively  different  possible  motions.  This  in  turn  means  that  we  cannot  develop  a  set  of  sim- 
plified composition  rules  as  we  did  for  fixed  axis  mechanisms. 

The  previous  observations  suggest  that  we  need  to  review  the  nature  of  the  qualitative  state- 
ments we  want  to  make  about  movable  axis  mechanisms.  There  are  several  possibilities  that  one 
can  consider.  One  is  to  find  an  interesting  class  of  mechanisms,  more  general  than  fixed  axis 
mechanisms,  for  which  it  is  possible  to  define  a  small  set  of  qualitatively  different  possible  paths 
that  cover  the  set  of  all  motions  and  that  are  well-behaved  with  respect  to  an  augmented  composi- 
tion operation.  Is  this  achievable  by  limiting  the  number  of  parts  in  a  mechanism  that  do  not  move 
along  or  around  fixed  axis?  A  second  approach  would  be  to  find  a  "coarser''  characterization  for 
qualitatively  different  motions  so  that  the  set  of  qualitatively  different  motions  becomes  finite. 

5.  Conclusions 

This  report  provides  a  formal  framework  for  reasoning  about  the  kinematic  behavior  of 
mechanisms  based  on  the  notion  of  configuration  spaces.  We  show  that  there  is  a  direct  relation- 
ship between  the  nature  of  the  configuration  space  and  the  possible  motions  that  objects  can  have. 
In  this  framework,  we  prove  that  the  configuration  space  of  a  mechanism  can  be  obtained  from 
the  composition  of  pairwise  configuration  spaces  that  correspond  to  the  kinematic  pairs  of  the 
mechanism.  We  prove  that,  to  analyze  a  mechanism,  it  is  sufficient  to  examine  all  pairs  of  objects 
that  can  be  in  contact,  compute  their  pairwise  configuration  space,  and  then  compose  them  to 
obtain  the  configuration  space  of  the  mechanism.  This  provides  a  justification  for  the  composition 
methods  presented  in  [Gelsey  87J  and  [Joskowicz  87]. 

For  a  special  class  of  mechanisms,  called  fixed  axis  mechanisms,  objects  can  have  only  five 
types  of  simple  motions.  Pairwise  configuration  spaces  can  thus  be  divided  into  regions  that 
correspond  to  different  types  of  possible  simple  motions.  We  show  that  the  composition  of  pair- 
wise  configuration  spaces  for  fixed  axis  mechanisms  can  be  carried  using  a  small  set  of  possible 
motion  intersection  rules.  A  complete  algorithm  that  determines  the  order  in  which  paiwise  confi- 
guration spaces  are  computed  and  composed  is  also  provided.  This  algorithm  is  based  on  the  con- 
cept of  constraint  propagation,  and  determines  the  order  and  the  direction  in  which  the  constraints 
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are  updated  and  propagated.  The  output  of  this  algorithm  is  a  region  diagram,  in  which  every 
region  corresponds  to  a  qualitatively  different  type  of  kinematic  behavior.  The  region  diagram 
contains  all  the  possible  kinematic  behaviors  of  the  mechanism,  and  thus  corresponds  to  the  result 
of  an  envisioning  of  the  mechanism.  The  region  diagram  is  similar  to  the  "place  vocabulary' 
presented  in  [Faltings  87]. 

The  region  diagram  can  be  used  to  determine  the  result  of  a  sequence  of  motions  applied  to 
different  parts  of  the  mechanism.  For  a  given  input  motion  sequence  we  show  how  to  construct  a 
state  diagram  containing  actual  motions  of  the  parts  in  the  mechanism,  the  different  kinematic 
behaviors  of  the  mechanism  and  the  transitions  between  these  behaviors. 

Finally,  we  indicate  that  the  simple  version  of  composition  operation  cannot  be  extended  to 
mechanisms  in  general  (movable  axis  mechanisms).  An  open  question  is  how  to  extend  the  class 
of  mechanisms  so  that  an  augmented  version  of  the  composition  operation  can  be  used  to  analyze 
them. 
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Figure  1.1: 


Three  different   geometries   for  a  prismatic  joint. 
The  configuration  space  of  the  slider  has  the  three 
orientational  parameters  constant  and  the  three 
positional  parameters  forming  a  line. 
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igure  1.2: 


Three  position  and  three  orientation  parameters  define 
the  exact  position  and  orientation  of  an  object  in 
three  dimensional  space 
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gure  1.3: 


Five  possible  simple  motions 
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Figure  2.1:      A  double  slider 
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Figure  2.2:      The  six  lower  pairs 


Figure  2.3: 


A  movable  axis  mechanism. 
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fure   2.4: 


The   Intersection  of  two  two -ri  <«„„«,., 
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Figure  2.5:      The  Clucth  System 
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Figure   2.6:  the    configuration    space    of   A    and    B    along    axis    O. 
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Figure  2.6: 
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A  Rod  and  Block  combination.  The  rod  can  only  rotate 
around  the  fixed  axis  01  and  the  block  can  only 
trans latealong  axis  02. 
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Figure  2.7:  Two  possible  partitions  of  the  Block  and  Rod 

configuration  space. 
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